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f^ , The reformulation of field theory for avoiding self-energy parts in 

■^ ' the dynamical evolution has been applied successfully in the framework 

^^ , of the Lee model, ||l| enabling a kinetic extension of the description. 

The basic ingredient is the recognition of these self-energy parts. [|| 
Mh, The original reversible description is embedded in the new one and 

appears now as a restricted class of initial conditions, g] This program 
H ^ is realized here in the reduced formalism for a scalar field, interacting 

r-{ ' with a two-level atom, beyond the usual rotating wave approximation. 

^^, The kinetic evolution operator, previously surmised, B is here derived 

from first principles, justifying the usual practice in optics where the 
common use of the so-called pole approximation!^ should no longer 
p\^ • be viewed as an approximation but as an alternative description in 

M , the appropriate formalism. That model illustrates how some dressing 

of the atomic levels (and vertices), through an appropriate operator, 
finds its place naturally into the new formalism since the bare and 
dressed ground states do no longer coincide. Moreover, finite velocity 
for field propagation is now possible in all cases, without the presence 
of precursors for multiple detections. 
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1 Introduction 

Quantum optics is not an ab initio theory and requires a quantum modeliza- 
tion of the interacting atoms and fields. The natural point of view starts 
with an atom described by its energy levels and a dipolar interaction with 
the field, in a time reversible formalism, using hermitian Hamiltonians.jP] 
Moreover, a welcome simplifying approximation (the so called rotating wave 
approximation) provides often an excellent approximation of the dynamics 
involved. Nevertheless, some problems subsists: they involve the descrip- 
tion of the instability of the excited levels and the acausal behaviour in 
energy transfer between atoms or precursors in the case of a double 
photodetection.p] The first problem is tackled by the use of irreversible 
elements inside a reversible quantum mechanical description. Indeed, for 
instance, the explicit attribution inside the Hamiltonian of a lifetime to an 
atomic excited level is the natural way to take into account its unstability. [^] 
The introduction of elements of a phenomenological origin into a microscopic 
description has then led to tremendous success in that field. The problems 
linked with the non-hermiticity are avoided by skilled use of the formalism. 
That practice has so far found no fundamental justification. 

In a previous paper, Q we have analysed in details a renowned paper |^] 
and we have shown that in fact. Hollow's approximations were equivalent to 
the use of kinetic equations to provide the description of the system. Implic- 
itly, such kinetic equations are naturally at the level of reduced distribution 
functions for the field. Indeed, when an arbitrary number of photons can 
be emitted from any initial state (except of course the true ground state), 
the most adequate description is a reduced formalism, applied currently for 
describing atoms in a fluid where quantal reduced distribution functions 
obey the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy Q in place of 
the Liouville-von Neumann equation. 

The use of those (kinetic-like) equations in ordinary quantum mechanics 
can be criticized on two bases. The first one is that they do not belong to the 
(reversible) framework of ordinary quantum mechanics and their theoretical 
justification is still missing (despite their extraordinary success). The sec- 
ond one is that they are always presented as the result of a largely justified 
approximation and not as an intrinsic property of the system. It seems un- 
satisfactory that a radical change of the formalism (a transition from a time 
reversal invariant to an irreversible one) results simply from approximations. 

Physicists associate the concept of an unstable state to an object that, 
outside external influence, would decay in a purely exponential way. Such 



an unstable state is beyond the reach of ordinary quantum mechanics and 
the tentatives to define it make trouble for normalization properties for 
instance. pH] The previous attempts for obtaining a purely exponential de- 



cay [11 1 were unsatisfactory since they did not incorporate the possibility 
of an excitation mechanism. The main aim of this paper is to show that 
an association between intuitive concepts and quantal description is possi- 
ble, through the single subdynamics approach, for the interaction of a field 
with a two-level atom, including the consideration of the counter-rotating 
terms, keeping in mind the generalization to multilevel systems. Acausal 
behaviours can moreover be excluded by an appropriate choice of a dress- 
ing operator. Under compatibility conditions to be satisfied at initial times, 
it is possible to show that the time reversal invariant and kinetic descrip- 
tions bring simply different aspects to light while keeping their equivalence 
(A spectral representation shares also that property. The derivation of a 
spectral-like representation for the Lee model through our single subdy- 
namics approach can be found in Ref. [p^]). 

A theory of subdynamics has been introduced thirty years ago by the 
Brussels group (see e.g. [13|, [^) for a dynamics provided by the Liouville- 



von Neuman equation. In that quest for the introduction of irreversibility 
inside the formalization of dynamics, the subdynamics concept has been 
shown to be fruitful. Different realizations are possible according to a choice 
of the vacuum, i. e. the choice of the degrees of freedom that are included in 
the resulting dynamics, while the other elements of the description become 
functional of the vacuum ones. A setback of that approach is thus a limita- 
tion on the class of possible initial conditions since they have to belong to 
the subdynamics. Therefore, the initial formulation should be in some sense 
overcomplete and contain degrees on freedom on which no control is possible. 
A way out has been the introduction of a transformation theory, intensively 



studied. 1 14] In the original ambition, an association of a so-called physical 
representation with real (energy-conserving) processes between renormalized 
(dressed) quantitities. It could not be carried out in a general and consis- 
tent fashion because of difficulties that have been reviewed, [p!^] Moreover, 
the problematics of ensuring positivity of the density matrix is lacking in 
late papers ||l^ where the positive character of the density matrices is no 
more listed among the requirements on the A transformation. 

In order to keep the completeness of the description, a new vacuum 
concept has been introduced |2[ for models in field theory, irrespectively of 
their classical [^] or quantum character. It is based on a dynamical anal- 
ysis of all possible contributions to the formal solution of the Liouville-von 



Neumann equations. To avoid the previous trap, the so-called single sub- 
dynamics approach Q is based on the existence of self-energy contributions 
to the dynamics. We can indeed accept that no control is possible on these 
processes. By definition, all self-energy parts have to be excluded from the 
vacuum. Their recognition implies that the initial dynamics is first extended 
by discriminating among the degrees of freedom according to their status 
with respect to preparation and observation. A subdynamics of the (ex- 
tended) dynamics is then introduced such that it encompasses the original 
dynamics but does no longer contain dressing processes. In that way, we 
obtain a reformulation of field theory that excludes self-energy contributions 
in the dynamics. They are now driven by the other degrees of freedom and 
provided by time independent functionals of the other degrees of freedom 
that are the motor of the evolution. Therefore, the same mathematical tool 
(subdynamics) as the Brussels-Austin group is used, but with a different 
realization, leading to different physical content, although a similar aim is 
pursued. 

In a previous paper in collaboration with C. George, y| we have been 
dealing with the Friedrichs model, equivalent to a specific sector of a two- 
level atom interacting with a scalar field within the rotating wave approxi- 
mation (RWA). The Friedrichs-Lee model has been treated in two different 
approaches. In the first one,|^] the existence of sectors has been used to 
perform the explicit construction of all the elements of the subdynamics 
super-operator inside the first non trivial sector. We have shown that a 
kinetic description exists that provides an exact and complete alternative 
to the (time reversal invariant) Schrodinger description. It is obtained by 
a double operation: An enlargement of dynamics, that enables the recogni- 
tion of self-energy parts, followed by the use of an appropriate subdynamics 
projector. The resolution of the model enables the explicit verification of all 
the claimed properties. That proof is welcome since the very existence of 
the concept of subdynamics has been questioned for instance by P. Coveney 
and O. Penrose. [|lq] Their argument is the incompatibility between branch 
points (generating "long time tails" ) and the kinetic description, preventing 
the subdynamics to provide the asymptotic behaviour. A general analy- 
sis of the situation, leading to their refutation, can be found in a previous 
paper. |15| In the second approach, ||l| the reduced formalism has been used 
to go beyond the limitation induced by the equivalence with the original de- 
scription. The focus has been on the subdynamics evolution super-operator 
and we have shown the possibility of extending quantum theory in a satisfac- 
tory way. Positivity and normalization are automatically ensured in the new 



kinetic description, starting from the reduced formalism. Such a construc- 
tion rests entirely on the existence of poles, independent of branch points, 
for diagonal matrix elements of the Green's function associated with the 
Hamiltonian. For this model, a clear-cut separation of poles and cuts does 
exist, leading to an intuitive description and justifying, through a derivation, 
a phenomenological approach, y] 

To show the robustness of that intuitive description, we treat here a 
more general system, consisting of a two-level atom in interaction with a 
scalar field, keeping the counter-rotating terms. The system under study is 
interesting from a physical point of view since the existence of long time tails 
(non exponential contributions) and deviations from exponential behaviours 
for very short times has long been recognized in it. On the other hand, 
physicists analysing experiments in optics are accustomed to use the so- 
called poles approximations. How can those empirical rules be justified 
from first principles? Are they valid only in some approximate way or do 
they fully reflect the physical reality. 

The associated subdynamics has to be constructed for the new system, 
outside the RWA. Such an extension involves formal modifications in the 
treatment, but the philosophy is the same. First of all, an analysis of the 
property of the Liouvillian using sectors [y] is no longer valid: A "reduced 
formalism" is therefore naturally required. All sectors are now coupled and 
have to be considered together. As a consequence, the possibility of a com- 
plete resolution of the model is lost with respect to the RWA case. 

Prom the expression of vacuum to vacuum elements of the resolvent 



of the generalized Liouvillian, 1 19] we have investigated whether a pole can 
be associated with all the matrix elements that do not involve neither an 
incoming nor an outgoing field particle, the generalization of the property 
for other kinds of matrix elements being straightforward. We have proven 



[19| that the notion of poles associated with the stable and unstable states 
is still relevant for the model under consideration treated inside the reduced 
formalism. 

The construction of the evolution and projection super-operator for a 
subdynamics rests on an analysis of the kind of behaviour for the contribu- 
tion of each term in a perturbation expansion of the resolvent (associated 
with the generator of motion) and the selection of the relevant behaviour. 

When compatibility (or equivalence) conditions are satisfied, [^ the ki- 
netic description provides somehow a change of representation, analog to 
a change of basis in standard quantum mechanics. Therefore, it will al- 
ways be possible to transfer the information on the system from the usual 



density matrix operator (or its analog for the reduced description) to the 
kinetic description and vice versa. Moreover, the intuitive way of looking at 
the system is recovered in the kinetic description in terms of incident field, 
outgoing field, dressed atomic levels (including their attribute of unstabil- 
ity represented by their lifetime, see Ref. |^ for instance). Therefore, the 
common practice in optics will be justified from first principle and will no 
longer be the result of approximations. Let us underline the analogy of the 
present approach with the theory of renormalization. Indeed, in both theo- 
ries, the aim is to take properly into account the self-energy contributions. 
The renormalization theory tackles the problem of removing ultraviolet di- 
vergences at the level of the wave function and to derive finite corrections. tl 
The description of an unstable state requires the introduction of an imag- 
inary part in the renormalized energy: The original framework cannot be 
respected. In the present theory, in terms of a reduced density operator, no 
such problem arises. 

In this effective realization of our approach on a less simple non trivial 
example, we will not bother about the formal properties of the subdynamics 
super-operator that have been well established [^ and do not depend on the 
particular realization of the choice of the vacuum. We will not worry either 
of the compatibility conditions: we know that they do exist Q and reflect 
the possible equivalence between the original and the kinetic descriptions. 
We focus on the derivation and the properties of the kinetic evolution super- 
operator, as defined by the single subdynamics approach. [^ The form of the 
obtained evolution generator leads us naturally to the introdution of dressing 



technique introduced long ago into the Brussels approach.] 14, 11| A dressing 
operator enables to shape the form of the evolution generator corresponding 
to an intuitive description, out of the reach of the original time reversal 
invariant description, that accepts only unitary transformations. 

In Ref. y], a reduced formalism has been proposed to treat the inter- 
action of a two- level atom with the electromagnetic field. Various kinetic 
equations have been accordingly surmised and justified on physical grounds. 
However, their derivation from first principles was outside the scope of that 
paper and we intend here to fulfill that missing part: In this way, we pursue 
the progressive introduction of the characteristics of the method to the ef- 
fective description of the interaction between the electromagnetic fields and 



^It has been recently noticed [[L7| that the apphcation of the s ing le subdynamics ap- 
proach to a well adapted formalism of classical electrodynamics bw leads directly to a 
divergence-free description, without any need to a substraction procedure. 



atoms. 

In §2, the reduced formalism for the system is briefly recalled and the ex- 
tension of the dynamics (distinction between the various photons) is treated. 
The formal properties of the subdynamics super-operator are briefly cited. 

§3 is devoted to the elements of the subdynamics super-operator that 
enable the obtention of the kinetic operator for the elements describing the 
population of the atomic levels. They do not involve physical photons The 
elements for the atomic dipolar moments are considered in §4. 

The main difference with respect to RWA is already apparent in the 
elements of the evolution generator 9 that do not involve photons (the purely 
atomic part of the evolution operator). Inside the RWA, the bare and dressed 
ground state coincide and the dressed excited state is then derived directly 
by the construction of the kinetic operator. Outside the RWA, the structure 
of the kinetic operator is no longer the same: both states that appear in it are 
susceptible of evolution. Indeed, already in usual quantum mechanics, we 
know that a change of basis is required with respect to the bare states outside 
RWA. Therefore, the stable ground state does not coincide with the state 
described by the kinetic operator. On physical grounds, it is required the 
true ground state being time independent and the excited state decaying. 
Since we are dealing in a reduced formalism, we have to translate such 
properties into similar ones for the matrix elements of the reduced density 
operator, in order to define a physical representation. 

Through a dressing operator, usual in the subdynamics approach, |jl^, ^] 
a procedure enables to fix the problem in §5 in a completely satisfactory way. 
After dressing, the structures of the kinetic operators inside and outside the 
RWA are common. Let us note that the atomic model under consideration 
is the first example treated for which a dressing operator is required. Our 
dressing, in the single subdynamics approach formulated in reduced formal- 
ism, is not equivalent of defining new states in an Hilbert space formulation, 
such as in [pO|]. 

For completeness, the (somehow trivial) effect of passive photons is ex- 
plicitly treated in §6. The one physical photon absorption process is also 
considered in the same §6 and leads, as in the case of the Lee model, O 
to the introduction of dressed interactions between the atom and the field. 
The one physical photon emission process does not present unexpected new 
features. 

The vertex dressing is considered in §7. We discuss the change in the 
interaction that can be induced by the dressing operator: Strict causality in 



the exchange of photons between atoms can be ensured in ah cases, without 
the usual presence of precursors due to a finite lower bound in the energy 
spectrum. [^-||2^ The use of a kinetic description, as opposed to a reversible 
one, is a main ingredient to allow that property, beyond the reach of an 
hermitian generator of motion. 

Concluding remarks are presented in the last part of this paper. 

2 The model - Reduction - Indiscernability - 
Extended Dynamics - Subdynamics 

The model considered is a two-level atomic system interacting with a field, 
without considering the rotating wave approximation nor explicitly the mo- 
mentum change due the recoil of the atom. The Hamiltonian of the system 
can be written as 

H = uJia\^ai + uJoa^aQ + '^ijOka^oik 

k 

+ H(^i|ofca5*"ao(afe + aA^) + V'ofc|ia(l"ai(afc + afc)) • (2.1) 



In another paper, [|19[ we have established the analytical properties of some 
Green's functions associated with the Hamiltonian ( | 2.1J ). Since the Hamil- 
tonian enables an arbitrary number of field particles being present in the 
future of any state, due to the counter-rotating terms, we are in a situa- 
tion similar to that met in statistical mechanics and a reduced formalism is 
required for the description of the degrees of freedom associated with the 
field. In the computation of the mean value of observables, such a formalism 
replaces the trace operation by a vacuum expectation value for the field. It 
has been developed in extenso in Ref. Q and recalled in Ref. j|l[ and we 
will be satisfied with a reminder of the main features without entering into 
details. 

The many body system of interest is described by a density operator p 
that obeys the Liouville-von Neuman equation 

idtp = Lp, ( 2.2) 

where the Liouvillian L is the commutator with the Hamiltonian H given by 



( |2.l|) . In terms of factorizable superoperators (AxB), defined as {AxB)p = 
ApB, it is given by 

L = H X I -I X H, ( 2.3) 



where / is the identity operator. A (factorizable) superoperator (A x B) will 
be called a connecting superoperator if both A and B are different from the 
identity operator. A superoperator such as in L that acts with an operator 
on one side of the density operator and with the identity operator on the 
other side will be called a non-connecting superoperator. It does not prevent 
from writing the formal solution of equation ( | 2.1| ) for p under a factorized 
form (exp —iHt) p exp iHt. 

The associated reduced density operator p obeys an equation 

idtp = Lp. ( 2.4) 

The unperturbed part Lq of L is the same as the unperturbed part Lq of L. 
The potential dependent part Ly of L contains, in addition to the potential 
dependent part Ly, new connecting contributions L'y given by 

L'y = Y^ (Vi\okafao + Vok\ia^aij x a^ 
k 

- X!"*: ^ (^ilofcO^^ao + V'ofciiaJ'ai) • ( 2.5) 

k 

If we note by a Roman letter the states of the Hilbert space (that letter 
defines the state of the atom and the wave numbers associated with the 
photons present) and if we use the notation < a\p\b >= pab, the evolution 
equations take the form: 

idtPab{t) = Y Lab.cdPcd{t), idtPab{t) = ^ Lab.cdpcd{t) ( 2.6) 

c,d c,d 

with an obvious definition for the matrix elements of the evolution super- 
operators L and L. 

A formal solution of these equations is provided with the help of inverse 
Laplace transform as 

Pc^bit) = E^/^^e"'*f7^) . fcdiO), (2.7) 

c,d 



2vri j-y Vz — L/jj^cd 



P^bit) = E^/^^e— *f^)^ p,M, (2. 

c,d 



2vrW7 \Z ~ L/ ab.cd 



^We keep the notations of Ref. |l| where the operators inside the reduction proce- 
dure bears a bar accent M while the symbol tilde is introduced to refer to the extended 
dynamics. With respect to Ref. pi, the reduced density operator is noted p in place of a. 



where the path 7 hes above the real axis. These forms enable immediately 
perturbation expansions in terms of the potential V. In the first part of this 
study, we limit ourselves to the case of the couples (ab) and (cd) referring 
to diagonal discrete states without field particles. The resolvent of L that 
plays a role in ( | 2.71 ) can be written as a convolution product of appropriate 
resolvents of H. 

Analytical properties of the resolvent of L and L in (| 2.7l - | 2.8| ) have been 
studied extensively in Ref. |jT^ (in particular the elements ((z — L)~^)^-^ ^-^, 

{{z - -^^)"^)oo.oo)' ((^ - ^)~^)ii.ii and {{z - i)"^)oo.oo)- ^he main techni- 
cal point involved was the ability of distinguishing the effects of the branch 
points and of the poles required for the construction of the single subdy- 
namics. In the reduced formalism, the analytic structure of the resolvent of 



the evolution super-operator can no longer be examined |19] a priori via a 
convolution involving the Green's functions associated with the Hamiltonian 
operator and this has far reaching consequences on the level of the analytic 
properties. 

We recall some of the properties that have been established. 



The Green's functions associated with the Hamiltonian ( 2.1 ) are noted 
as in the RWA case 

r]i{z) \z-HJii rjo{z) \z-HJqq 

Inside the RWA, the r]o{z) function reduces to {z — uq). "Bar" functions 77 
represent Green's functions associated with the operator {—H) 



1 \ 1/1 



rji{z) \z + H J 11' rio{z) \z + H 



( 2.10) 
00 



so that the relation ( | 2.7| ), for a = b = c = d=loia = b = c = d = 
involves the resolvent Rii.ii{z) and Roo.oo{z) given by 

-1 /■ , 1 1 

-Kii.iU^;) = —^dU' 



2-Ki J^i r]i{u) r]i{z — u)^ 
Roomiz) = ^ fdu ) . _ / . . (2.11) 



The r){z), r]{z) functions are required in ( 2.11] ) for a positive imaginary part 



of their argument (0 < ^7 < Qz). Therefore, those functions will be analyt- 
ically contined from above into the lower half plane Qz < 0. When branch 
points are met, the cuts will be placed parallelly to the imaginary axis. All 



functions of complex arguments that we shall introduce share that property 
and we shall dispense them from a upper index "+" that would indicate the 
way the analytical continuations are performed. As a consequence of that 
convention, for instance, f/i (z) for '^z < cannot be computed directly from 
( 2.10D but requires the computation of ((z + H)~^)-^-^ for 3=2: > and then 



an analytical continuation. It is usual to note that fact by writing fii{z) 
but to avoid too cumbersome notations, we use an implicit convention to 
keep the notation fii{z). The same convention holds for other functions to 
be introduced soon. 

The T]~^ functions are analytic in the upper halfplane Qz > and we 
first analyse their singularities in the lower halfplane 9z < 0. The properties 
of the fj functions are a translation of those of the r] functions. From the 



convolution product ( 2.11 ), the singularities of the resolvent into the lower 
halfplane Q'z < are then inferred. 

From a perturbation expansion with respect to the interaction V, the 
rji function ( | 2.9| ) can be written in general as ryf = J2'^=oii^ ~ -f^o)^^)ii 
(/i(z)((z — i/o)"^)!!)" where fi{z) represents the sum of all the contribu- 
tions leading from state "1" to itself with a condition of irreductibility with 
respect to "1" (all intermediate states implied in the sum have to be ^ of 

"1") : 

oo 

/l(^) = E ((^((^ - Hor'WrrlV)ll. ( 2.12) 

The index irrl recalls the restriction on the intermediate states. The diag- 
onal matrix elements of {z — Hq)~^ are called propagators. 

A visual representation of the contributions can be obtained using a 
diagrammatic representation, such as the one developped in appendix B of 
Ref. 1^] for the same system. It corresponds in the present case to the 
Feynman diagrams. We have then trivially 

11 11 , , 

( 2.13) 



rji{z) z-uJi- fi{z)' rjo{z) z-lvq- fo{z)' 

The summations present in the expression of the functions / have to be 
understood in the infinite volume limit where, for instance, J^k l^ilofcP . . . — > 
/q°° dLOv'^{uj) . . . We assume that the function v'^{lo) vanishes as uj for small 
values of its argument (This may be understood as it is composed of a factor 
k'^ arising from the jacobian to obtain spherical coordinates and the square 
of a usual factor -j=^ arising from the matrix element of the potential.). 



10 



The thesis estabhshed in Ref. [19| is the fohowing: Under this behaviour 
of the potential matrix elements Vi|ofc, V"o|ifc, Vo|ifc, Vi|ofc for small value of 
the argument k , we have the consistency of the following statements for the 
singularities of the functions r] and /. 

l)The function (r/i(z))~^ presents a pole at some point z = toi + C and a 
well defined residue ^i at that point. 

2)The function {r]Q(z))^^ presents a pole at some point z = loq + 6 and a 
well defined residue ^o at that point. 

3)The functions fi{z) and /o(-z) present logarithmic singularities at the 
points z = LOi + C and z = loq + 5. 

4)The function fi{z) behaves as {z — loi — C)^ ln(2; — toi — (^) in the vicinity 
of its singular at point z = uji + ( and as (z — ujq — 6) ln(z — ujq — 6) in the 
vicinity of its singular point z = loq + 5. 

5)The function fo{z) behaves as (z — ujq — 5)^ ln(2; — uiq — 5) in the vicinity 
of its singular point z = ujq + 5 and as (z — uji — Q ln(z — uJi — Q in the 
vicinity of its singular point z = loi + (. 

(^ has an negative imaginary part while 6 is real. Inside the RWA Q, the pole 
of {r]i{z))~^ has been noted 6i. Here, we bring to the fore the displacement 
with respect to the unperturbed value coi. For the "bar" functions, we have 
the similar properties: the poles of [fjiiz) — 1) and {fJQ{z))~^ are resp. at 
z = —ioi + C and z = — wq + S, with 6 = —6, iC, = {K)* ■ 

The explicit demonstration has required a self-consistent analysis of the 
various contributions. For the system we consider here (outside the RWA), 
the previous proof of the existence of the poles inside the RWA [^ has 
been adapted: the structure of the Hamiltonian Green's functions is no 
longer as simple and the diagonal elements of the Green's functions present 
simultaneously a pole and a branch point for a same value in the complex 
plane z, while poles and branchpoints do not coincide inside the RWA. 

An analysis of that new intertwining has enabled to desantangle it in 
order to be able to make a statement about the existence of poles indepen- 
dently of a choice of a particular Riemann sheet. The problematics to which 
an answer has been given is thus the existence of a "pole approximation" 
§ beyond the RWA. 

From the convolution form, ( | 2.11| ) it is then established that Ru.iii^) 
and -Roo.oo(-z) have respectively a pole aX z = Q + C, and z = Q. In a similar 
way, it can be easily established that Rw.iq{z) and Rqix)i(z) have simple 
poles respectively iov z = uji + Q — ujq + 5 and z = ujq + 5 — uji + Q. 

That analysis could not be reproduced mutatis mutandis. Indeed, the 
new generator of motion in the reduced formalism is no longer likewise sim- 
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ply connected with the Haniiltonian since the notion of reduction leads out- 
side the haniiltonian formalism. Therefore, the existence of "connecting ver- 
tices" ffl prevents a direct analysis in terms of a convolution involving matrix 
elements of the Green's function associated with the Haniiltonian. Neverthe- 
less, a further analysis has shown that both matrix elements -Rii.ii(z) and 
^00.00 (-2) of the Green's function associated with the reduced Liouvillian L 
have poles z = Q + C, and 2 = 0, with well defined residues, jl^ That result is 
not unxpected on physical ground: the natural time behaviour of the atom 
should not be modified by the reduction procedure. The existence of those 
poles is a necessary requisite for the construction of the subdynamics. 

The indiscernability of the field and its consequences on the computa- 
tions has already be treated in Ref. [l| and it has been noted that they are 
minimal. We will not repeat the same considerations. Indeed, in Ref. [Q, 
care has been taken on that question. 

The evolution for field and atom reduced distribution functions is gov- 
erned by an evolution operator L. Following the approach initiated with 
the potential scattering [^, ^] and the atom in interaction with the field 
inside RWA,|||] an extended dynamics is now introduced, that rests on a 
distinction between all degrees of freedom, once the self-energy parts have 
been recognized. We use the same notations as in Ref. [||] , a Roman letter 
I to represent a photon involved in the self-energy contribution, a Greek 
letter ^ for a photon leaving the atom (emitted photon) and a Greek letter 
A for an incident photon. The number of reduced distribution functions to 
be considered is multiplied accordingly and the evolution is then governed 
by an evolution operator L: 

L = ujio^ai X I — I X LOia'^ai + ujQa^ao x I — I x woo^ao 

+ ^ujialai X I - I x^ujialai 
I I 

A X fJ- /i 



+ 



+ 



Y [^m^i"'o{ai + a-t) + VoiliaQaiiai + alYj x I 
I 

I >^Y1 [^m^to-oi^^i + o-l) + Vo/|ia^ai(a/ + a^ 

Y (^i|OAa]^aoaA + Vox\ia^aiax) x / 
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A 

+ Z! (^i|0/.Oi'aoa+ + VQ^\ia^aia+^ x / 

+ J2 {^mo-tao + Vo/|ia+ai) x a+ -J^f^i^^ (^i|o/ai'«o + Vo/|iat}'ai) 
+ Yl {^iloxatao + VoAiiflo'ai) x a^ - ^ oa x (Viioxatao + Vox\ia^ai 

A A 

( 2.14) 

The system is now described by a set of new reduced distribution functions 
p related to the extended dynamics. The constitutive relations connects the 
original set p to the new one p.[^ Indeed, the new description contains obvi- 
ously more degrees of freedom than the original one. Since the observables 
are defined originally in terms of p, we have to specify the relation between 
the two descriptions. The constitutive relations do not involve the "/" pho- 
tons (involved in a self-energy process) and enable the interference between 
emitted and incident photons. 

The construction of the subdynamics requires a classification of the states 
into two classes, the vacuum states and the correlated states. They are 



obtained traditionnally |14] by the introduction of an superoperator P that 
projects on the so called vacuum states. In our case, correlated states contain 
at least one intermediate field line (of the / type). Vacuum states contain 
thus only incoming and outgoing field lines. 

The construction rule for the subdynamics operator can be formulated 
on the formal solution of the evolution equation in the extended dynamics: 

-1 r . _.,. / 1 



Pabit) = T.^J/'<^'"' (^) PcdiO). ( 2.15) 



The path c is chosen above the real axis for t > 0. Analytic continuation of 
the integrand from Qz > to Qz < is therefore required upon integration 
over z. The symmetry with respect to time inversion is thus broken by the 
procedure. All integrations on intermediate field lines are (at least) formally 
performed to permit an analytical continuation from above of the functions 
of z so defined, placing all cuts parallely to the imaginary axis. The rule is to 
pick up the contribution of the poles associated with vacuum states. The just 

13 



mentioned analytical continuation from above enables to avoid accidental 
coincidence of the poles associated with the vacuum and correlation states. 
The physical poles are located at a value defined by the poles associated 
with the atomic states, i.e. z = ( + CjZ = 0, z = uji + C — coq + S, z = 
loq + 5 — uJi + Cj ill addition to frequencies associated with physical field 
lines (incoming or outgoing). For instance, the vacuum- vacuum element 

( {z — L)^^ ) may have the following physical poles: z = (^ + ^ — lox, 

V / IXOfi.lXl 

Z = -UJx, Z = OJi + C - LOo + 5 - UJx + OJfj,, Z = UJo + 5 - LOi + C - UJx + UJfi- 

In the perturbation approach, we have the following expansion for the 
resolvent Rah.cd{z) of L (involving the resolvent R^{z) of the unperturbed 
hamiltonian Hq and the interaction part Ly of the liouvillian) : 



RabUz) = E (^°(^) \LvR\z)T) ( 2.16) 

•^ — ' \ L J / ab.cd 

n=0 

For the resolvent Rab.cd{z) of L we have similarly the expansion: 

oo 

Rab.cdiz) = E (^°(^) \LvR\z)T) ^ . ( 2.17) 

•^ — ' \ L J / ab.cd 

n=0 

The singularities of the resolvent R{z) of L are defined on resummed expres- 
sions. Therefore, useful expressions are obtained by considering irreductible 
operators with respect to the vacuum. For the study of Rab.cdiz), where 
no extension of dynamics has been defined, we may define here the set of 
vacuum states by the states without any field particles .H Such a vacuum is 
useful since it will provide a point of comparison for the elements of Rab.cd 
that do involve neither incident nor emitted photons. The collision operator 
[14 1 il'iz), is defined as the sum of irreductible fragments leading from a 



vacuum state to another: all intermediate states have to imply at least one 
field particle: 

oo 

^l^ab.cdiz) = E {\LvR\z)TLv) ^ ^^. ^ ( 2.18) 

'^ — ' \L J / ab.cdiirr] 



n=0 



where the index irr means that al intermediate states involve at least one 
field line. 



^In the more usual approach by the Brussels group, the vacuum has been defined by 
the set of diagonal elements or by an adequate extension in the approach by the patterns 
of correlation. M 
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We express also the perturbation expansion of the resolvent R{z) in 
terms of irreductible operators: For the complete Liouvillian L, we introduce 
by analogy the irreductible operators Wab.cd (that could also be noted by 

i'ab.cdiz)) 



Wab.cd{z) 



oo 

E 

n=0 



LvR{z) 



Lv ] ( 2.19) 

ab.cd{irr) 



When a, 5, c, d represents states without photons, the operators Rab.cd{z) 
and Rab.cd{z) coincide. In order to be able to make the connection between 
the singularities of R{z) and R{z), let us define irreductible operators that 
involve at least one of the connecting vertices L'y ( 2.5) as: 



Tab.cdiz) = (\LvR''{z)r Lv) ^ ^^. ^ (2.20) 

V L J / ab.cd(trr,con) 

where the subscript con implies that at least one of the vertices L is a 
connecting contribution L'y ( 2.5| ). The link with the above introduced 



operators V' and W is then obvious: 

Wab.cd{z) = 1p{z)ab.cd + T{z)ab.cd ( 2.21) 

In view of the characteristics of the vacuum states (an incoming field line 
cannot be reintroduced and an outgoing field line never disappears), the 
contributions can be easily classified according to the number of vacuum 
field lines involved and computed in a recurrent way. Therefore, we first 
consider the contribution of terms that do not involve any field line. 

3 Elements of E without field 

Our aim in this section is the computation of the elements of subdynamics 
super-operator Tjab.cd{t) when a, 6, c, d represent states without photons. 
As in the case of the RWA,[j|] they determine the evolution super-operator 
Qat that does not involve absorption nor emission processes. These elements 
of Tiab.cd{t) are computed from the corresponding elements of the resolvent. 
As no physical field lines are involved, the extension of dynamics plays no 
role: only virtual photons, involved in self-energy contributions, play a role: 
We may use the Liouvillian L in place of L. The elements of the operators 
Rab.cd{z) and Rab.cd{z) Coincide and the analytical properties of Rab.cd{z) 
have been established in Ref. ||l^ for some diagonal-diagonal elements and 
the extension of the proof to the other matrix elements is straightforward. 
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The elements of S are expressed in terms of the irreductible operators 
W, ■0, T. For the elements of these operators without photons, some useful 
relations may be considered. The last vertex (at the extreme left) in the 
perturbative expansion ( 2.19) for Wab.cdiz) corresponds to the absorption 



of the last virtual photon. For that vertex, the replacement of the element 
of L'y by the element of Ly that absorbs also the last photon, or vice versa, 
transforms the contribution to Wab.cdiz) into a contribution to W with other 
values of the first two atomic indices, leading to the relations: 

Wn.cd + Woo,,d = 0, (3.1) 

Wio.cd + Woi.cd = 0. (3.2) 

Since ^11.00(2;), tpoo.ii{z), ^10.01(2) and ^01.10(2;) vanish, the links Q 2.211 ), 
with the previously introduced operators tp and T, provides the following 
relations 

Tiim{z) = -^oo.oo(^) = -V'oo.oo(z) -7bo.oo(^;), 

700.11(2) = -Wn,u{z) = -iJn.n{z)-Tu.ii(.z), (3.3) 

^10. 01 = —Woi,Ol{z) = — V'01.01 — ^01. 01; 

Toi.io = -^^10.10(2) = -V'lo.io -T'lo.io- (3.4) 

This section is devoted to the elements Sn.n, Sqo.oO) ^11.00 j ^ioo.ii of the 
subdynamics super-operator T,(t) in order to obtain the elements Gn.n, 
©00.00) ©ii.ooi ©00.11 of the evolution operator G. The other elements with- 
out physical field lines, such as Sio.iO) are not dynamically connected to the 
previous and are treated in the next §4. 

The relevant elements of S are the same as those of Sn.n, Sqo.oO) ^ii.oOj 
^00.11 computed by choosing the diagonal atomic states as the only elements 
of the vacuum. In terms of the irreductible operators, the following relations 
hold: 

Ru.uiz) = .Rll. 11 + -Rll.llTii.iii?!!. 11 + i?ii. iiTii.oo.RoO.il) 

Roo.uiz) = .Roo.ooTbo.ii.Rii.li + -Roo.oo^oo.oo-Roo.ii) 
-Rii.oo(2) = .Rii. 11^^11. 11.R11. 00 + -Rii.iiT'ii.oo-Roo.oO) 
-Roo.oo(2) = .Roo.oo + -Roo.oo^oo.ii-Rii. 00 + ^00.00^00.00^00.00- (3.5) 

For instance, the first one expresses the possibility of transition from (11) 
to itself either directly (without connecting vertices), or with at least one 
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connecting vertex that involves Tn.n or Tuqq. The system (| 3.5) can be 



solved easily, leading to [cf. Ref. [|19[, §4] 



-Rii.ii(^;) 
-Roo.oo(^) 
-Rii.oo(^) 
Roo.u{z) 



z 



Woo.ooiz) 



z{z- Wu,uiz) - VFoo.oo(z)) ' 

z-Wn.ii{z) 
z{z-Wu.iiiz)-Woomiz))' 

Tumjz) 

^(^ -H^ii.ii(^)- ^00.00 (^))' 
100.11(2;) 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



z{z - Wii.iiiz) - Woomiz)) ' 
The relevant analytical properties of those expressions have been studied in 



Ref. |19|. The coincidence of the pole z = 9 with the corresponding pole at 
z = 6 for the Liouvillian L has been established. The residues at those poles 
are well defined, although the functions ^^11.11(2;) and Woo.oo(-2) present also 
branch points for the values of z = and z = 9 = 9 = C + C,- 



'^aa.bbit) IS then computed from 



■'aa.bb 



(t) 



-1 

27ri 



dz e-'^'R 



aa.bb 



iz) 



-1 

27d 



dze 



-izt 



L 



( 3.10) 



aa.bb 



where the path contains only the two relevant physical poles: The first one 
is the obvious z = 0, the other one is a pole at z = 9{= 9 = C + 0- The 
value 9 satisfies formally the equation 



9 = Wu.u{9) + Woo.ooi9), 

the solution does not depend of the chosen Riemann sheet. 
We have by direct integration 



( 3.11) 



Sii.ii(t) 



W, 



00.00 



(0) 



(VFii.ii(O) + 1^00.00(0)) 



exp —i9t \Ai\ 



,W^ 



11.11 



{0) 



9 



3.12) 



The residue at the pole z = 9 has been noted as |^ip. That manner of 
writing has been introduced by analogy with the similar computation inside 
the rotating wave approximation ^ and does not imply a possible existence 
of ^1 for instance. Indeed, we have shown in Ref. |19| the existence of \Ai\'^ 
as a sum of residues and not as a product of separate contributions. 

The main difference with the treatment inside the RWA can be seen 
on that expression: two poles play a role in the expression of this (doubly) 
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diagonal-diagonal element of S. Therefore, we have no longer the association 
of a temporal behaviour with a state. Inside the RWA, the bare ground 
state coincides with the true ground state while we have now to establish 
the connection. 

Let us make some comments on the relative values of the contributions 
in a perturbation expansion in powers of the potential V. Indeed, it is easily 
realized that VFoo.oo(O) vanishes at the second order while VFii.ii(O) is finite 
at that order and coincides with the second order value of 9. This difference 
of behaviour has its origin in the stability of the state at the lowest order. 
Therefore, the second term starts as 1 in a perturbation expansion in the 
potential while the first one is of order V^. 

The other elements of S are evaluated in appendix A. The vacuum- 
vacuum elements of the super-operator S(t) for the value t = define the 
operator usually called A in the subdynamics theory. Its derivation and the 
obtention of its inverse is now straighforward. The corresponding elements 
of the evolution operator in the vacuum subdynamics can be obtained 
from the relation (P is the projector on the vacuum elements): 

pt{t)P = exp-ietA. (3.13) 

Therefore, by multiplying the time derivative of S(t) at t = by the inverse 
of the operator A, we get: 



eii.ii 

©00.00 
©00.11 



Wu.uiO) 



{Wn.u{0) + Woomm 
i ^00.00 (0) 

(VFii.ii(0) + Vroo.oo(0))' 
g H^ii.ii(O) 

(W^ii.ii(0) + Woo.oo(0))' 



A - o ^00.00 (0) /-QiA^ 

°"-°° - "'(^ii.ii(o) + Tyoo.oo(o))- ^'-''^ 



Some obvious comments are a direct consequence of the form ( 3.14 ) of the 
@ operator. A first remarkable point is that this result does not obviouly 
depend on a choice of a Riemann sheet: that expression requires only the 
existence of the pole 9: the z = values for the elements of the irredictible 
operator W are obviously well defined. That operator G preserves the norm 
of p, thanks the obvious relations ©n.ii = — ©oo.ii and ©n.oo = — ©oo.oo- 
One of the eigenvalue of © vanishes while the other one is the sum 9 of its 
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diagonal elements. A unitary transformation can lead to a Jordan form, 
with the cancelation of the elements $11.00 ^md $00.00 of a new evolution 
operator $. The structure obtained inside the rotating wave approximation 
would then be rederived in an exact way. This point will be considered (see 
§5) after the examination of the properties of the other elements of S that 
do not involve field particles and describe atomic dipolar moments 

4 The atomic dipolar moment in the reduced for- 
malism 

We are interested in this section in the elements without field Sio.iO) ^^oi.oij 
Sio.oi, Jioi.io that are not dynamically connected with the previous ones. 
Those elements play a role in determining the properties of the dipolar 
moment of the atom. As in the previous section, the extension of dynamics 
plays no role, the index "tilde" can be dropped and express the elements in 
the reduced formalism in terms of the elements in the unreduced one's. Our 
starting expressions are therefore the formal solutions [formulae (| 2.7 - 2l8|) ] 



for the evolution of p and p, where the couples {ah) and {cd) refer to the 
off diagonal discrete states (10) and (01). This can be easily symbolized by 
using a notation [aa): a is one when a is and vice versa. An element of 
comparison is provided by studying first R{z). That resolvent of L that plays 
a role in ( |2.7| ) can be written as a convolution of the corresponding resolvents 
of H and —H. Since the interaction involves the change of occupation 
numbers of field particles by one unit for each absorption or emission process, 
the states a and c on one hand, h and d on the other hand have to be 
identical. Only diagonal elements of the resolvent of H and —H are thus to 
be considered and we have 

Ra-a.a-a{z) = {^) =^f duf^—) ( ^-— 

\Z-LJaa.aa 27rz Jo<S«<Sz \u - H J ^a \Z - U + H y ^a 

(4.1) 

involving the functions r/ and rj given in ( | 2.9[]~2T0| ). The analysis of the 



convolution enables [24, 1£] to ascertain the presence of an isolated pole to 
both diagonal elements of R. In terms of the corresponding poles of the 
Green's functions associated with the Hamiltonians H and —H (see §2 or 
Ref. ||l9|), they are given respectively hy z = ui + C — ^0 + S with a residue 
^1^0 for (ad) = (10) and hy z = loq + 5 — toi + C with a residue ^o-^i foi" 
the other case (ad) = (01). 
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The computation of the elements of S is detailed in appendix C. The 
poles that play a role are still located at z = ui + C — loq + 6 = 6io and 
z = ujQ + 5 — ui + ( = 5oi. They now appear together as poles of each 
element. The result can be symbolized as: 

y _ _|'+'l — p~*'5lO*n' - --Lp~^^Olta _ _ 
'-'aa.aay^ ) — '^ '-"■aa.aa \ ^ yaa.aaj 

^aa.aayt ) — 6 0!aa.aa ~r C Paa.aai \ 4./J 

where the values of a and /3 can be read on the equations ( | C.Sl) . In a 
perturbation expansion, aio.io and /?oi.oi start as 1 while the other elements 
behave at least as V"^, since they involve connecting vertices. 

The corresponding elements of the evolution operator Q in the vacuum 
subdynamics can also be obtained from the relation ( | 3.13| ) through the 
previous procedure (|^d| is defined in ( C.12[ )) 



^aa.aa — -l. yOwOlaa.aa ~r OQ\Paa.aa)\P^aa.aa ~r Paa.aa)i 

Qaa.aa = {Sioaaa.aa + SoiPaa.aa){aaa.aa + Paa.aa)- (4.3) 

In opposition with the case met in the RWA, the operator O is not diagonal 
inside the sectors defined by the atomic dipolar moments. 

The value of I e~ ) _ _ can be obtained directly from the knowledge 

V J aa.bb 



of the elements of S and A : 




\ / aa.bb 


= 


(M ') (t) 

\ / aa.bb 






= 


(e-''^°*aaa.io + e-'^°^*/3aa.io) 


/ ^-1 

A 
\ 




+ 


{e-'^'''aa-a.oi + e-''^°^*/3aa.oi) 


A 
\ 



10.66 



01.66 



(4.4) 



We know from the expression of $](t) that the two eigenvalues of that op- 
erator are b\Q and 5oi for the elements under consideration. Therefore, 
an invertible transformation can provide a diagonal evolution operator that 
has these two values as the only non- vanishing elements. 
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5 Physical representation for the purely atomic 
part of p 

Denoting Pat{t) the elements of p that do not involve photons, in the absence 

- at 

of an incident field, we have, by definition of Q , the set of equations: 

d ~ -at^ 

-p^,{t) = -ie p,t{t), (5.1) 

at 



where the relevant elements of the operators are given in ( 3.14| ) for 



the diagonal-diagonal elements and in (| 4.3) for the off-diagonal ones. Let 



us examine first the diagonal-diagonal elements. The elements 0oo.oo and 
©ii.oo would vanish if Woo.oo(O) = 0. We may check explicitly this prop- 
erty in a perturbation expansion. The second order value W^qq oo(0) involves 
only contributions that are present in the unreduced formalism and van- 
ishes accordingly. We have therefore to consider the fourth order value 
Wqqqq(O). After some algebra, a non- vanishing contribution, that involves 
the connecting vertices L'y ( | 2.5| ) is obtained (see Appendix B for details of 
computation) . 

^oaoo(O) = -27riY.Y.\^iiok\'\Viiok' 

k k' 



^SiuJi - 


-UJQ - 




1 


"''^-1- 


-ujQ+ujk'y' 








(5.2) 



Since we have 

Wi?,n{0) = -27riJ2\Vi\0k\^S{ui - oj^ - Uk) = 0^^\ ( 5.3) 

k 

the first non-vanishing contribution to Ooo.oo is of the fourth order 

^oaoo = 0^'^ E I^Ho^'I't ^T S2- ( 5-4) 

The denominator is positive defined and no regularisation ("ie") appears. 

= (4) 

Obviously, Oqo.oo does not vanish. As a consequence, poo does not describe 
the atom in its ground state. For comparison, Qn n starts at the second 

= (2) - - 

order and we have On.n = 0^'^' . We know from the expressions of S that the 
operator 0, relevant for the diagonal-diagonal elements has two eigenvalues, 
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and 0. If we take into account the conservation of the norm, the operator 
O can be related to an operator $ that has the structure met inside the 
RWA: 

$11.11 = 0, $00.11 = -6, $00.00 = $11.00 = 0. ( 5.5) 

Let us note that the natural structure for the evolution operator does not 
require an hermitian operator (or star-hermitian), in opposition to the early 



attempts in the subdynamics approach.] 14] We define a new reduced den- 
sity operator p (in the so-called physical representation) connected to the 
original one p through an invertible dressing operator x 

f{t) = x-^m, m = xf{t), (5.6) 

and such that the evolution of p is governed by the operator ^. Therefore, 
the following relations have to hold: 

<^ = X-'§X, § = x$X"'- (5.7) 

~ —P 

If we impose that the traces of p and p are the same, we have moreover 

(x"^)ii.ii + (x"^)oo.ii = 1, (x""^)oo.oo + (x""^)ii.oo = 1, 

Xii.ii + Xoo.ii = 1) Xoo.oo + Xii.oo = 1- (5.8) 



We have therefore to determine x by imposing conditions d 5.8| ) and d 5.7 ) 



The last conditions (5.7) can be made explicit, taking (| 5.5|) into account: 



©aa.fefe = Xaa.ll$ll. 11 (X ^)ll.feb + Xaa.00$00.1l(X ^)ll.bb- (5.9) 

The values of the elements of x ^re computed from these conditions in 
Appendix D and the unambiguous result is 



Xii.ii = Xoo.oo 

Xii.oo = Xoo.ii 

(x^ )ii.ii = (x^ )oo.oo 



Wii.ii(O) 



VFii.ii(0) + Woo.oo(0)' 

M^oo.oo(O) 
VFii.ii(0) + ^oo.oo(0)' 

1^11.11(0) 
W^ii.ii(0)-Woo.oo(0)' 



/ -1\ r -1\ W^OO.Oo(O) / c 1^^ 

(X )oo.ii = (x )ii.oo = -w,,,m-Woo.oo{Oy ^'-''^ 
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The form of the operator ^ and the condition on the norm determines 
thus entirely the dressing operator x for these elements. The two diagonal 
elements of p can thus be interpreted as the probability of finding the atom 
resp. in the ground and excited levels (cf. |jl|). They are normalized and 
have the expected associated "free" motion (in the absence of an incident 
field). 

We now turn to the elements describing the dipolar atomic momentum. 
We know from the expressions of S that the operator Q has two complex 
eigenvalues inside the considered subspace, namely Siq and (Jqi- The operator 
G should be related to the operator $, the elements of which are given by: 



$ 



10.10 



no, 



$ 



01.01 



501, 



$ 



01.10 



$ 



10.01 



0. 



( 5.11) 



The natural structure for these elements evolution operator requires here an 

hermitian operator, as opposed to the case of the diagonal elements. The 

—P 
relevant elements of the reduced density operator p are connected with the 

original one p through the same relation ( | 5.6| ) and its evolution is governed 

by the operator <1> ( | 5.11 ). Therefore, the same relations ( 5.7 ) have to hold 

inside the sector, but no condition on the trace appears. 

Denoting by \xd\ the determinant of the x matrix in the dipolar sectors, 

we have from the inversion of a standard 2x2 matrix 



vX Jaa.aa — Xaa.o 



\Xd\ 



{X-')c 



-Xa 



\Xd\ 



( 5.12) 



We have therefore to determine x by imposing conditions (5.7) that can be 



made explicit, taking ( 5.11 ) into account 



b 



Qaa.cc — 7 ,XanM"bh\X 



5.13) 



We know ( 4.3 ) the explicit value of the elements of @, but we will not 
proceed by direct identification of ( 5.13| ) and ( 4.3| ). In order to obtain the 
values af the elements of x, it seems easier to compare the expressions of 

'^ and from its expression (| 4. 



exp —iQt deduced from x(exp —i^t)x~^ and from its expression ( | 4.4| ). We 
proceed in appendix D with the identification of the terms with the same 
exponential behaviour {exp—i5iot or exp—i6oit). The elements of x can 
be determined up to a free parameter. A remaining indetermination is not 



new inside the context of subdynamics.|14| It does not affect the evolution 
equations but concerns the relation between p and p. That indetermination 
is linked to a choice of relative phase between the basic states used for the 
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description. The dressing operator is indeed apt to change the original 
choice. We can fix it by connecting, for instance, the value of |x| with that 
of \A£)\ in a "usual way": |xP = \-^d\ or introduce a criterion that ensures 
the positivity of the density operator if required. We will not elaborate here 
further on this point. 

For the elements considered, the reduced density operator p satisfies 
the kinetic equations that we have postulated in Ref. [^. That derivation 
is exact (no approximation has been required). The kinetic equations ap- 
pear naturally when the physical representation is introduced. As already 
stressed in Ref. [Q], physicists in optics use them naturally, without realiz- 
ing their profound origin, while considering conceptually such equations as 
useful approximated equations arising from some pole approximation. The 
formulation of observables in the formalism p is more natural that in the 
original representation. What we call the ground state and excited states of 
the atom are the objects that behave as the diagonal elements of p . in the 
absence of an incident field, the excited state decays in a purely exponential 
way while the evolution of the ground state is only due to the transfer arising 
from the decay of the excited state. Such a decay takes place with the life- 
time as it is usually computed from 5-matrix or Green's function formalism. 
Our formalism, that eliminates the reabsorption of the field by its source, 
leads automatically to a kinetic description while preserving the equivalence 
with the original Liouville-von Neumann equation, providing the compati- 
bility conditions are satisfied at the initial time.pl The atomic observables, 
computable with pat (= Pat by the constitutive relations), have also to be 
modified by the dressing operator (x or more precisely x~^) if they are de- 
fined in the original representation, so that the conservation of their mean 
value is ensured. If the observables are defined by the kinetic properties, 
such as the probability of finding the atom in the ground or excited state, 
they can be expressed directly in the physical representation. Our approach 
enables moreover a consistent departure from the original time-reversible 
description while preserving normalizability and positivity. [Q] The element 
Pii has a clear physical meaning as the probability of finding the atom in the 
excited state, irrrespectively of the state of the field. Inside RWA, in |j|], we 
have shown that an initial condition describing the atom in the excited state 
(only pii does not vanish) does not belong to an admissible initial condition 
if the equivalence conditions are satisfied. Such an initial condition is how- 
ever fit to discuss normalizability and positivity in the irreversible departure 
from quantum mechanics but its relevance in physics has to be established 
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on an experimental basis. 

6 One incident field 

6.1 One passive incident field 

We look first for the description of the elements of the evolution operators 
Q and $ involving an incident field line that is not absorbed. They can be 
obtained by considering the elements of T,{t) that involve one incident field 
line at the extreme right and the same field line at the left. The atomic 
states are explicit in our notations and are represented by the letters a, b,c, 
d (a, 6, c, d represent the atomic state or 1), while the present photons 
are written at the right of the atomic state. We will consider therefore the 
following elements of S: SaAfe.cAd; ^abx.cdx- The couple of atomic states {ab) 
and (cd) involved at the right [cd) and at the left {ah) of the elements of 
S are of the same nature to get a non-vanishing contribution, namely a 
diagonal couple or an off diagonal couple on both sides. 

Our starting point is the formal expression of the contributions to pa\b{t) 
arising from the part of initial conditions involving only one incident field 
line: 

-Paxm - E ^ / ^^ ^""' {-^^ /^^Ad(0) ( 6.1) 

c,d -'" \z-LJ^xb.c\d 

The detailed calculation can be found in Appendix E and provides the ex- 
pected result. If we define lab.cd as 1 when a = c and b = d {lab.cd vanishes 
for the other possibilities) we obtain 

©aAfe.cAd = ^\Iab.cd + ^ab.cd ^abX.cdX = —^xlab.cd + ^ab.cd ( 6.2) 

The dressing does not depend on the incident photon and the natural choice 
for the dressing operator x is 

XaXb.cXd = Xab.cd XabX.cdX = Xab.cd ( 6.3) 

and therefore 

^aXb.cXd = ^xlab.cd + ^ab.cd ^abX.cdX = —^xhb.cd + ^ab.cd ( 6.4) 

The generalization of that property in the case of the presence of an arbitrary 
numbers of field lines (incident or emitted) is obvious as long as they do not 
interact with the atomic variables. 
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6.2 One absorbed incident field line 

We look for the elements of the evolution operators Q and ^ that describe 
the absorption of one incident field line. They can be obtained by considering 
the elements of S(i) that involve one incident field at the extreme right and 
no field at the left. We consider therefore the following elements of S (a, b, 
c, d represent the atomic states or 1) Tiah.cXd, ^ab.cdx- The couple of atomic 
states (ab) and (cd) involved at the right (cd) and at the left (ab) of S are of 
different natures: We have a diagonal couple on one side and an off-diagonal 
couple on the other side. The intrinsic evolution of that kind of couples 
have been studied in previous sections. If the diagonal couple is at the 
left, we have to consider: Sn.iAO, Sqo.oai, Sh.oai, Sqo.iao, Sh.iqa, Sqo.oia, 
S^ii.oiA, ^ioo.ioA- Similar elements have to be considered in the case where 
the diagonal elements are at the right: Sio.iAi, Siq.oao, Sqi.iai, Sqi.oao, 

Slo.llA) Sio.oOA; ^Ol.llA) Sqi.ooA- 

Our starting point is the formal expression of the contributions to Pab{t) 
arising from the part of initial conditions involving only one incident field 



-Pab{t) - E^/^- 



^~izt 



c,d,>- 

.(6.5) 



PcXd{0)+( ^1 Pcdx{0) 



^^ ~ L/ ab.cXd \Z~L/ab.cdX 

The details of the computation can be found in Appendix E and takes 

the form: 

~ ~ ~- ^-1 ^ ^-1 ~ ~ ~-i 

^aa.bXb = 2^^^ ''^'^■^'^^^ jcc.dd^dd.eXeK^ Jee.bS + Z^(©^)aa-cAc(^ Icc.bb 
c4,e c 



-1 ,^.?~-, .~—l. 

Jee.bb + Z- 
d,e c 



2_^&aa.ddAdd.eXe{A )^gi,i, + }_^{eA)aa.cXc{A )cg.bfe. (6.6) 



Similar expressions hold for the other elements @aa.bbXi ^aa.bXb-, Qaa.bbx- 
The derivation of the elements describing an emission process can be 

performed in a completely similar way. We will not dwell on these terms 

since they do not introduce new ideas or properties. 

These expressions may serve to determine the corresponding elements of 

the dressing operator x and of the evolution operator $ (see next §7). 
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7 Vertex dressing 

We are interested in the possibilities for the evolution operator $ that de- 
scribes the absorption of one incident field line. We use our previous knowl- 
edge of the elements of Xi that do not involve the field. These elements have 
been determined in previous sections, using physical requirements. For the 
computation of the element of <1> describing the disparition of the incident 



field line, we may use the link ( 5.7) between Q and $ and introduce a new 
element Xgh.cXd for the dressing operator to get 



ab.cXd 

= Z^ (X )ab.efQef.ghXgh.cXd + 2^ (X )ab.ef&ef.gXhXgh.cXd 
ej,9,h ej,g,h 

+ 2^ (?<: )ab.eXfQeXf.gXhXgXh.cXd- ( 7.1) 

Let us call <!>"* the part of <1> that does not involve transitions in the field. 
The elements of <1>'^* have been determined in §5 [ (| 5.5| ) and ( 5.11 )] and 



in §6 d 6.4^ ). Using this known value of $"*, the property XgXh.cXd = Xgh.cd 
(see ( 1 6. 3D ) and the value of the element {x~^)ab.eXf of the inverse of x, that 
can be computed in the same manner as the corresponding element of the 
inverse of A in appendix E, we get 

^ab.cXd = 2_/ (^" )ab.ef&ef.gXhXgh.cd+ 2_/ ^ab.efiX~ )ef.ghXgh.cXd 
ej,g,h ej,g,h 

~ Yl iX~^)ab.efXef.gXh^fxh.cXd- ( 7-2) 

ej,g,h 

Defining the X operator such that its only non-vanishing elements are 

Xab.cXd = 2^(X )ab.efXef.cXd, 
e,f 

^ab.cdX = Z^{X~ )ab.efXef.cdX, ( 7.3) 



we can express the last two terms in ( 7.2) as a commutator of X with the 
operator ^^^. 

^ab.cXd = 2^ iX~ )ab.ef&ef.gXhXgh.cd+ {[^'^ ,X])ab.cXd 
ej,g,h 

^ab.cdX = Z^ iX )ab.efQef.ghXXgh.cd+ {[^'^ ,X])ab.cdX- (7.4) 
ej,g,h 
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The elements ^ab.cXb and ^ab.cdx are expressed as the corresponding elements 

of 0, with a dressing bearing on the atomic levels, plus the commutator of 

$"* with an undetermined X operator. The role of that indetermination can 

be understood from the introduction of the dressed reduced density operator 
-P 



p (I 5.6). Indeed, the determination of the elements of X is equivalent 



to that of X [from ( 7.3)]. The presence of a non-vanishing X operator 



^p 
means that the elements of the physical reduced density operator p^^ are 

also dressed by the elements p^^b and Pabx- The dressing of the atomic 

states involves the incident photons. This is by no means mandatory: we 

can choose to dress the states by the self-field only. Therefore, the future 

determination of the X operator requires a personal choice of the basic 

states for describing the atom.[P Such a choice cannot be implied by the 

formalism. 

Due to the structure of the commutator, the contribution to $ depending 

on X vanishes for the resonnance processes when the width of the states is 

neglected: This contribution plays a role in the off resonnance processes. 

This point can be illustrated by considering some first order non-vanishing 

contributions to $, the element 'J'^x OAi ^"-"^ instance. We can replace in 



( 7.4 ) the X and x ^ operator by unity. The first order contribution to 
©ii.OAi is nothing but T4|oa- '^ii.ii is at least from the second order {O^"^') 
while ^oAi.oAi provides at the lowest order a contribution independent of 
the coupling (^q^^^ qai = "^o + i^a — "^i)- Therefore, 

^ILOAl = ^l|OA - (^0 + ^A - a;i)Xii.0Al. ( 7.5) 

If Xii.oAi is not singular for lox = uji — ujq, its value is irrelevant for deter- 
mining 'I'll oAi ^^ ^^^ ^^'^^ energy resonnance {ujx = uji — loq). 

It has been long noticed |25|-[^] that the kind of coupling (| 2.1]) be- 



tween atoms and field is not completely satisfactory with respect to the 
causal propagation of the field, i.e. its finite propagation. Indeed, in the 
problem of transfer of excitation from one atom to another, if the measure- 
ment concerns also the state of the initially excited atom (in a case of non 
local observables, as treated in case I of Ref. [^) precursors appear and 
causality is not strictly respected. [^ The same kind of acausal behaviour is 
met Q in the description of a double photodetection of the light from an 
atom admitting a cascade decayQ. Those non causal behaviours are linked 



^In Ref. P9|, causality is claimed to be respected in the case of a double photodetection 
but that property is limited to two-level systems while multilevel systems were considered 
in Ref. § 
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to the finite lower bound in the spectrum of the exchanged photon and dis- 
appear when the integration over the energy of the photons is extended to 
— oo. This was common practice. Various authors pOj-pSl have claimed 
to have succeeded in proving causality of propagation of light : They all 
used at some stage an equivalent procedure. In the case of |33|, for in- 



stance, the sign of a contribution is changed on the base of its smallness.[]g] 
In later works, that include the counter-rotating terms, that procedure has 
no longer been required for proving strict causality in either the framework 



of two-level systems for all observables,[34, 35,|22] or, for multilevel systems 
for local observables.[^ Nevertheless, the general case, beyond the two-level 
systems and local observables, still contains acausal behaviours, even when 
the counter-rotating contributions are considered.]^, |^ 

In Ref. p, a suggestion has been proposed to ensure a strict causality in 
all possible cases by proposing new terms of interaction between the atom 
and the field. The present formalism is more appropriate for realizing that 
possibility of including in a simple way causality into the equations of mo- 
tion without modifying the energy spectrum of the field. With respect to 
the usual approach, we have indeed a supplementary degree of freedom since 
the symmetry beween absorption and emission is lost.O, |l| In the study of 
the coupling of the quantum field with a quasiclassical source, Q it has been 
established that an appropriate association of super-operators for the cre- 
ation of the field \£c in Ref. Q] and absorption process \£ in Ref. 0]] leads 
to a strictly causal propagation, without precursor : The retarded solution 
for the field propagation appears automatically, irrespectively of the atomic 
state. We can use the possibility of dressing to meet the causality require- 
ments. As we have seen, the dressing bears on non-resonant contributions. 
Therefore, we can require that the $ elements, describing the interaction 
with an atom at some point r, have the appropriate structure for a local 
coupling (k is the wave number associated with A or ^) 

1 1 

ff, _ Moi>^ ^ik.r ?f. _ labs ^-ik.r 

^ab.cXd — (Pab.cd.—^^^ > ^ab.cdX — (Pab.cd—^^ 



if, A^Tn „— jk.r /f, Mf^m „«k.r ( ^ a\ 

^atib.cd = (Pab.cd^^e , "Sabt^.cd = -(Pab.cd^^e . ( 7.6) 



The dressing by x cannot change the value of for the resonant contribu- 
tion at the lowest order in the coupling but can provide any chosen value 
for the other (off-resonant) elements. It is thus perfectly possible to choose 
X such that these conditions are fulfilled. The change of relative sign be- 
ween the two tems describing absorption and emission is capital. When an 



29 



emission, by a pointlike atom at some point r, is combined with the absorp- 
tion by another one, at some other point ri, two processes lead to the same 



change in the atomic occupation numbers ( 7.6). The field emitted by an 
atom through ^afib.cd, and absorbed by the other atom through ^a'b'.c'Xd', 
provides a time dependence as exp(— iw^t) while the other process {^ab^i.cd 
and ^a'b'.c'd'x) involves exp(iu;fci). Both contributions can be combined and 
their sum is equivalent of having a domain of integration over the photon 
spectrum from — oo to -|-cx), but this procedure is now an integral part of the 
evolution equation. Therefore, the resulting contribution ignores the pres- 
ence of a finite lower bound in the spectrum of the field and causality can 
be fully respected, without precursors or other oddities:]^] The propagator 
corresponding to a retarded solution for an electromagnetic process appears 
naturally, corresponding to Ritz's point of view on the origin of irreversibil- 
ity, rather than Einstein's conception. [^] Of course, such a radical change 
in the form of the coupling has to take place while preserving positivity of 
the retrieved density operator. 

The emphasis on the physical interpretation of the elements of the de- 



scription is not the only difference with the Brussels- Austin group. |1€ 



Indeed, as can be seen on the Friedrichs model, |37] their approach is based on 
the consideration of Gamow vectors obtained through a generalized eigenval- 
ues problem, outside Hilbert space, providing complex values. That obten- 
tion have been inspired by the early attempts of constructing a sub dynamics. [11 
Besides problems linked to normalization (such states are of null norm and a 
set of left and right bicomplete and biorthogonal eigenvectors for the Hamil- 
tonian has to be introduced), and positivity (that is not a requirement of 
their A transformation), such states can only decay (by construction) and 
the description of an excitation of the atom by an incident wave packet is 
outside the possibility of that approachE In contrast, the presence, in our 



kinetic operator $, of non-diagonal elements in the field (see ( 7.6| ) for in- 
stance) enables a coupling between an incident wave packet and the atom. 
Therefore, the work of that Brussels- Austin group is more appropriate for 
an abtract discussion about irreversibility, p^] through the introduction of 
"diagonal singularities", than for the description of atom-field interactions 
required in quantum optics. 



■''That criticism can also be addressed to our previous work |ll|] and has been one of 
our motivations to depart from the early attemps of the Brussels group. 
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8 Concluding remarks 

This paper enables to understand the link between reversible and irreversible 
formulations of interacting atoms and fields. This later formulation requires 
a description in a Schrodinger type description. The existence to an in- 
vertible transformation enabling a transition from one description to the 
other one may be perceived as a surprise: when compatibility conditions are 
satisfied, both formulations are equivalent. The physical representation en- 
ables physicists to provide in quantum optics a description conform to their 
intuition, without any loss of generality or the introduction of approxima- 
tions. An atomic level is characterized by the values of its (dressed) energy 
and its lifetime. Dressing contributions are now excluded through in the 
structure of the evolution operator that involves the often used distinction 
between external and emitted photons. Quasiclassical state for the external 
field can be considered. |Q] The external field description naturally factorizes 
inside the reduced density operator. Acausal behaviours can now be ex- 
cluded through the choice of the dressing operator. A modelization for an 
ideal photodetection device P finds naturally his place in that framework. 
The relation with the original description involves the dressing operator, the 
constitutive relations and the compatibility conditions. A direct formula- 
tion of an initial condition for a problem inside the physical representation 
is nevertheless possible, with the possibility of dealing with an extension of 
quantum mechanics, as has been considered in Ref. [l|]. An initial condition 
where only the physical excited state is present cannot be obtained from 
an initial condition in the original formulation, [y, ^ The necessity of such 
an extension of quantum mechanics cannot be excluded a priori but should 
result from experiments. 

This paper has focused on the role of dressing inside the single subdy- 
namics approach. In the first papers within that approach, no need for a 
dressing has appeared. The main difference of the system inside and outside 
the RWA is indeed the following: In presence of counter-rotating terms, the 
obtention of the ground and excited states of the atom is no longer auto- 
matic and requires a dressing operator. Indeed, inside the RWA, the bare 
ground state can be identified with the true ground state but that property 
is lost outside the RWA: It will not convert itself into an invariant state 
and an adequate combination is required that enables the identification of 
the ground state and of the excited state through their temporal depen- 
dences. Renormalization through a dressing operator takes naturally place 
for dealing with physical ground and excited states, physical atomic dipolar 
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moments, pointlike interactions and causality. The recourse to these entities 
to make physical predictions depends of course of the observables consid- 
ered. If they are defined in the original representation, no need to go into 
the physical representation is present, except for the simplicity and trans- 
parency of the evolution in that representation. If the observables are to be 
precized, they are more naturally defined inside the new representation. In- 
deed, the precise meanings of the ground state, the excited state, the atomic 
dipolar moments resort to p where temporal behaviours enable the identi- 
fication. The recourse to the new description can also be justified according 
to the manner in which the initial conditions have to be formulated. If we 
intend to consider the scattering of a wave packet impinging on an atom in 
its ground state, we need to know how to describe correctly the atomic state 
and the present approach answers that question. 

We have not considered in this paper the constitutive and compatibility 
conditions since they appear here in a fashion similar to our previous papers. 
We are aware that positivity should require a more detailed analysis to place 
restrictions on the undetermined elements of the dressing operator. 

This paper shows morever how the modelization in optics is justified 
from first principles without the need of approximations such as the pole 
approximation. In our kinetic equations, the explicit attribution of a lifetime 
to an atomic excited level is the consequence of the original Liouville-von 
Neumann equations when going into an "historical" representation. 

The kinetic equations obtained in the present paper are at the level of 
reduced distribution functions for the field, and have been introduced in [Q|. 
The kinetic equations that were surmised in that paper are now properly 
derived. 

The dressing procedure introduced in the early approach of the subdy- 
namics theory has proved to be still relevant and fruitful in the present 
work for more elaborate systems, closer to quantum field theory. We are 
nevertheless still facing some basic indetermination corresponding to the ar- 
bitrariness in a choice of basis vectors. That indetermination is not met for 
the purely atomic part, except for a (trivial) choice of the relative phase 
between the ground and excited atomic states: The atomic dipolar moment 
is not completely fixed. The indetermination can modify the form of the 
evolution operator when absorption and emission processes are considered 
and can be used to ensure a strict causality, in spite of a finite lower bound 
in the energy spectrum. 

We have not dealt with simultaneous processes Q that are to be present 
both inside and outside the RWA, the simultaneous absorption and emission 
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process appears at a lower order outside the RWA than inside. Since they 
involve necessarily nonresonant contributions, they could be of the same 
importance as higher order resonant contributions, not included into the 
model Hamiltonian. 
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A Other elements of S between diagonal atomic 
states 

We now turn to the second element Sqo.oo- From the following relations 



obtained in Ref. |1£] 



^ii.ii(^;) + Room{z) = - H ttz ^^ ttj ^—r, ( Al) 

z z - Wu,u[z) - Woo.oo[z) 

we have 

f;ii.ii(t) + Soo.ooCi) = 1 + exp -m \Ai\\ ( A2) 

and therefore 

The first term starts as 1 in a perturbation expansion in the potential while 
the second one is of order V"^ . 

We now turn to the next element Sn.oo; computed from ( 3.8). The con- 



tributions to Sii.oo(i) emerge also from the poles at z = and z = 9. The 
existence of those poles, established in Ref. [l^l for Ilii.ii(t) and Iloo.oo(0 
enables to prove their existence for Sii.oo(t). We can easily see that the 
residues of those poles are well defined, in particular that the numerator 
is uniquely defined for that value. Indeed, by comparing thir perturbation 
expansion (the relevant part of the unpertubed propagator in reduced for- 
malism EP{z) coincides with EP{z))^ we have 



Rmm{z) = E (^°(^) [LvR\z) 



n=0 



00.00 
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n=0 



11.00 



(A4) 



If in the second expression (of Rii,oo{z)), we modify the last vertex we get 
a contribution to .Roo.ool-z) with a change of sign. The term n = 0, without 
vertex Ly, of .Roo.ool-^) cannot be recovered in that way and we have to add 
its contribution. Similar considerations hold als when the last two indices 
ar "11" in place of "00" and we have 



-^11.00(2;) 
-Roo.ii(^) 



■-Roo.oo(^) + - 
--^11.11(2;) + -, 



(A5) 
(A6) 



We can check the compatibility of the two expressions ( | 3.8 ) and ( | A5| ). 
Using the expression (| 3.7| ) for Roqoq{z) we get: 



R 



11.00 



[z) 



-Room{z) + 

-Woo.oo{z) 



W^ 



11.11 



(z) 



+ 



1 



z{z- Wii,ii{z) - ^00.00(2:)) 



Wu,u{z) - Woo,oo{z)) z 

(A7) 



If we take into account the relation ( | 3.4|) between Woq,oo{z) and Tiioq{z), 
we get the equivalence between the two forms ( | 3.8|) a nd ( | A5| ) for Riiqq{z). 
Its interest lies in the following remark: in Ref. |jl^, only the existence of 
the poles for the elements .Rii.ii(2:) and .Roo.oo(-z) has been considered. The 
formulae (| A5| ), ( A6 ) show directly that the poles of .Rii.oo(-2) and .Roo.ii(-z) 
are automatically well defined. We get therefore: 



5^11.00(0 - 
and similarly: 

Soo.ii(0 = 



-Tii.oo(O) 



(VFii.ii(0) + VFoo.oo(0)) 

-roo.ii(O) 
(VFii.ii(0) + 1^00.00(0)) 



+ exp —i9t \Ai 



2^ii.oo(^) 



(A8) 



+ exp-iet\Ai\'^TQo,ii{e). ( A9) 



An alternative expression for Sii.oo(i) and Soo.ii(i) can be obtained from ( 
"A5| ), ( [AeD or from ([TsD 



m.oo 



it) 



t^oo.oo(O) 



(W^ii.ii(0) + T^oo.oo(0)) 



exp —iOt \Ai 



:Woomi0) 



AlO) 
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f, M ^11.11 (0) -z^,, J ,2 ^11.11 W , ,-n, 

Let us note the obvious relations, that can be derived from ( | A5| ) and ( | A6| ) 

too,u{t) + tu.u{t) = l, too,oo{t) + tum{t) = l. ( A12) 

Those relations ( |3.12| ), ( [A3|) , (| AlOD and ( | All]) can be symbolized as: 



^aa.bb{t) = Oiaa.bb + e * Paa.bb, ( Al3) 

where the values of a and /? can be read on the preceding equations. In a 
perturbation expansion, an.n and /3oo.oo starts as 1 while the other elements 
behaves at least as V'^. 

The vacuum- vacuum elements of the super-operator S(t) for the value 
t = define the operator usually called A in the subdynamics theory. The 
corresponding elements are 



Aaa.bb{t) — Oiaa.bb + Paa.bb, (^14) 



or more explicitly 



T ^ ^OO.OO(O) I T ,2 ^11.11 W 

"■" (w^ii.ii(o) + T^oo.oo(o)) ' " e ' 

7 _ Vt^ll.ll(O) ^ I T 12 ^00.00 W 
^00.00 — TTTT /r,\ , jxr 77^ + -^1 ^ 

T ^ m)O.Oo(O) _ I T 1 2 ^00.00 W 

"•°° (W^ii.ii(0) + T^oo.oo(0)) ' " ^" ' 

I _ ^11.11 (0) |,-|2 ^11.1lW .,.,. 

The corresponding elements of the inverse A of the A operator can already 
be computed, independently of the elements of the subdynamics superop- 
erator involving the field by the inversion of a two by two matrix. The 
determinant Ad of that matrix is 

■Ad = ^11.11^00.00 — ^00.11^11.00 = l-^il , ( A16) 

as can be shown by direct computation using (| 3.1l| ). We have therefore 
directly 

X-'\ ^ ^ll.ll(O) I J 1-2 ^oo.oo(^) 

/ 11.11 (H^ii.ii(0) + VFoo.oo(0))' " 9 ' 
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/oo.oo Wi.ii(0) + P^oo.oo(0))' " ^ ^ ' 

rl\ ^ V^^OO.OO(O) , T ,-2 ^00.00 W 

Ai.oo (W^ii.ii(0) + Woo.oo(0))' " 0" ' 

/oo.ii {Wn.ii{0) + WoomiO)y " ^ '^ ^ 

B The fourth order contribution to Woo.oo(O) 



.(4) 



We derive in this appendix the expression ( | 5.2 ) for Wqqqq{0). We have to 



consider ah irreductible contributions leading from the diagonal state 00 to 
itself. We have 16 contributions to evaluate. Our convention to denote a 
matrix element < a\p\b >= pab of the (reduced) density operator p is to 
write for a and b in the first place the state of the atom. This convention 
avoids the need of a separator when writing the indices ah. Using that 
convention, the possible succession of correlated states to be considered to 



evaluate Wqo.oo from ( C.IO ) can be described in the following way. 



We first consider contributions that do not involve the connecting vertex 
L'y and are included in V'oo.oo- They correspond to the following successions 
that involve as last vertex (extreme left) a change of the first index: 
(00;lkO;Okk'0;lkO;00), (00;lkO;Okk'0;lk'0;00), (00;lkO;lklk';lkO;00), 
(00;lkO;lklk';01k';00). The corresponding contributions are called Ci{z), 
C2{z), 6*3(2), Ci{z). They are: 



1 1 

, i^^iiofcri^^iiofc'r : 

k k' 



^iW = EEl^nofcl'l^nofc'l' 



Z — UJl + UJQ — OJk Z — UJk — UJk' 



1 



Z — LOl + LOO — UJk 

C2{z) = EEl^Hofcl'I^Hofc'l 

k k' 

1 



2 1 



Z — LOl + LOq — ^k' 

Csiz) = EEl^Hofcl'l^Hofc'l 

k k' 

1 



Z — LOl + LOq — LOk Z — LOf^ — ^k' 



2 1 1 

Z — OOl + LOo — OOk Z — CVk + ^k' 



Z — LOl + LOq — LOk' 
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C^iz) = EEl^Hofcl'I^Hofc'l' 



1 



- UJl + UJo - UJk Z - iOk + iOk' 

1 



Z + UJl — LOq + UJk' 

(B.l) 

The contributions involving as last vertex a change of the second index 
(such as (00;01k;00kk';01k;00)) can be obtained form these expressions by 
replacing in the propagators uj -^ —uj for all values of the indices. The 
corresponding contributions are called Cy, C2', C3/, C41. The limit z ^ 
has to be considered to get the contributions to Wqq qq(0). For the terms Ci 
and C2, the limit z ^ can be taken in a harmless way since no propagator 
can be resonant. Therefore, we have 

cm = EEi^iiofci'i^iiofc'i' — -^ ^ -^ , 

kk^ -Wi +UJo-UJk -UJk - UJk' -UJl +UJo-UJk 

C2io) = EEl^i|o;^l'l^i|o;^'l^ T^ ^ T^ • 

kk^ -Wl +UJo-UJk -UJk - UJk' -UJl +UJO- UJk' 

(B.2) 

For obvious reasons, we have therefore: 

Ci(0) + Cr(0) = C2(0) + C2'(0) = ( B.3) 

We now turn to the contributions C3 and C4. The limit z ^ has to be 
taken carefully for the propagator (z — uJk + ^k')^^ since it is resonant for 
UJk = ojk'- Let us consider the sum C34 of C^ and C4 and take the limit 
z ^ whenever it is harmless: 

1 1 



C34(0) = limV5]|Vl|ofc|'|Vl|ofc''' 



""^^ k k^ Z- ^k+^k' -^1 +^Q-^k 

1 1 

+ 

-Wl + UJQ — UJk UJl— UJQ+ UJk' 

V Y^Y^IT/ |2|T/ |2 ^k' — ^k 1 

kk^ Z-^k+^k' -^1 +UjQ-UJk 

1 1 



—UJl + UJQ — UJk UJl — UJQ + UJk' 

EEi^iiofci i^iiofc'' 



UJl+ UJQ- UJk —UJl + UJQ - UJk UJl — UJQ + UJk' 

(B.4) 



37 



The limit z ^ has become harmless: the resonant factor plays no role 
since the fraction has a well defined limit. Therefore, we have anew the 
cancellation of (734(0) with 6*3/4/(0). 

We now turn to the contributions that involve the connecting vertex 
L'y. They correspond to the following successions that involve as last vertex 
(extreme left) a change of the second index: 

(00;Okl;Okk'0;lkO;00), (00;0kl;0kk'0;lk'0;00), (00;Okl;lklk';lkO;00), 
(00;0kl;lklk';01k';00). The corresponding contributions are called C^{z), 
Ceiz), Criz), Cs{z). They are: 

C5{z) = J2J2\^m\'\yi\0k'\'— — ^ ^ — 

V V Z + UJl - UJQ - UJk Z - UJk - UJk' 

1 

Oq — UJk 



Cj{z) 



Csiz) 



Z - UJl + LOo ■ 

>:>:\vm\ 
k k' 
1 


- Wfc' 

2|T/ 
1 yi\Ok' 


Z - UJl + UJo ■ 

k k' 
1 


- Wfc/' 

2|T/ 
1 l^llOfc' 


Z - UJl + LOo ■ 

>:>:\viiok\ 

k k' 
1 


- Wfc' 

2|T/ 
1 l^llOfc' 



Z + U!i — UJQ — UJk Z — UJk — UJk' 



1 



Z + UJl — UJQ - UJk Z — UJk + UJk' 



1 1 



Z + UJl - UJQ - UJk Z - UJk + UJk' 
Z + UJl — UJQ + UJk' 



B.5) 



The contributions involving as last vertex a change of the second index can 
be obtained form those expression by replacing in the propagators uj ^ —uj 
for all values of the indices. The corresponding contributions are called C5/, 
Cg/, C7/, Cg/. With respect to Ci{z), C2{z), 6*3(2:), 64(2;) , we note that in 
6*5(2), 6*6(2), 67(2), Cs{z) the first propagator is resonant. Let us consider 
the sum 656(2) of 65(2) and Cq{z) at the limit 2^0. Anew, the limit is 
taken whenever it is harmless. 



C56(0) = limV^ 1^110^1' ll^nofc'" 



1 1 



^"'^ kk^ Z + UJl - UJQ - UJk -UJk - UJk' 
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X ( i + = ). (B.6) 

\-UJi+UJo-UJk -Wi+Wo-Wfc'/ 

If we add the contribution C^'Q'iO), we get: 

""^^ k k' 

1 1 \ 1 

+ 



Z + LOi — UJo — LOk Z — OJi+ UJQ + UJk/ —^k — ^k' 

1 1 

+ 



^^^51 51 l^llOfcPl^llOfc'P'^K - ^0 - ^k) 

kk^ -^k-^k' 

1 1 

+ 



-UJl +UQ-UJk -^1 +^0- ^k' . 



k k' 
1 f 1 



— CJi + UJQ - UJk' \ 2(-^i + Wo) -1^1 + 1^0 - ^k' , 

(B.7) 

We have used the Dirac delta function to simphfy some propagators to get 
the last relation. 

We now turn to the sum Ci%{z) of Cj{z) and Cs{z) at the limit z ^ 0. 
Anew, the limit is taken whenever it is harmless. 



C78(0) = hmV^|Fi|ofeP|Fi|ofe,'2 



1 1 



^^° k ¥ Z + UJl - UJQ - OJk Z - UJk + i^k' 

1 1 

+ 

-LUi + LOq — LUk UJl — UJo + UJk' , 

kk^ Z + ^1 - ^0 - ^k Z - UJk + UJk' 

1 1 



-UJl +UJO- UJk UJl- UJo + UJk' 



1 

Z + UJli- 

1 1 



!i°^H5Il^i|0fel^l^i|0fe''^ 



° fc fc^ Z + tJl - CJo - tJfc 



-Wl + UJo — UJk UJl — UJQ + UJk' 



(B.8) 
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If we add the contribution C7'8'(0), we get 



C78(0) + C7'8'(0) = limVElFiiofel'lV^i 



|2 

„,^^^^ , .|.n,, , .|Ofc'| 

2— >U 

k k' 



1 1 

+ 



Z + UJi - UJo - UJk Z - UJl + UJQ + UJk 

1 1 

-UJl + UJQ - Uk UJl - UJQ + ^k' 

k k' 
1 1 

-Wi + (Jo — i^fc 1^1 — 1^0 + Wfc' 

fc fe' 
1 1 



2(— wi + LJo) wi — Wo + Wfc/ 

(B.9) 

If we combine the non-vanishing contributions, we are left with 

^oo.oo(o) = -27TiY,J2\^m\>m'\^^(^^-'^o-u^k) 



k k' 

1 



{UJI -UJQ +UJk')'' 



( B.IO) 



that is the looked-after expression. It is therefore manifest that VFqo.oo('-*) 
does not vanish. 



C Elements of S between off diagonal-ofF diago- 
nal atomic states 

The matrix elements of the resolvent of R can also be expressed in terms 
of the irreductible collision operator tp: Only elements ipaa.aa play a role 
while elements ipaa.aa and ipaa.aa vanish. We have the alternative forms with 



respect to ( | 4.1|) : 



Ra-a.a-aiz) = 7 ; T^- ( C.l) 

Z - UJa + UJa - iPad.ad{z) 
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The analytical properties of i^aa.aaiz) are determined by comparing the two 
expressions of Raa.aa{z) and are directly connected with the analytical prop- 
erties of the resolvents of H and —H. 

For the reduced formalism, the role of ■;/' is fulfilled by the operator W 
( 2.191 ), the two operators differing by T (\ 2.20| ) . The elements of R{z) are 



expressed in terms of these of R{z) and T{z) 

-t^aa.aa\Z) — t^aa.aa \ ^aa.aa-'-aa.aa^aa.aa i i^aa.aa^aa.aa^aa.aa^ 
^aa.aayZ) — t^aa.aaJ-aa.aa^aa.aa y ^aa.aa-'-aa.aa^aa.aa) \ ^•■^j 

A little algebra enables to solve this system of equations, using the relation 
( [2:21!) for Wab.cdiz): 






Z- UJa+UJa- Waa. 



aa 



{z-UJa+^^a- Waa.aa){z -Ua + UJa- Waa.aa) " Taa.aaTa^ 

T- - 

^aa.aa 



{z — UJa + UJa — Waa.aa)\Z — U}a + ^a — Waa.aa) — Taa.aaTaa.aa 

(C.3) 



This form shows that all elements share poles due to the common denomina- 
tor (we are interested only to these poles that go either to toi — wq or loq —uji 
as y^ — > 0, and not to the singular points arising from the numerators of ( 
C.3| )). The two off-diagonal elements {Raa.aa for a = 1 and a = 0) are also 
linked with the diagonal ones by (cf. ( | A5| )) 

{z - UJa + (^a - Ipaa.aa) Raa.aa = -{z - UJa + UJg. - Ipaa.aa) Raa.aa + 1 ( C.4) 

These relations can be interpreted by considering the last apparition, in 
Raa.aa , of an element of L'y, leading to a state without field particle (an 
odd number of L'y vertex has to be present), and its replacement by an 
element of Ly, leading to an even number of vertex L'y . The new expression 
can be placed in relation with contributions to Raa.aa- A supplementary 
contribution of Raa.aa has to be treated separately: Those relations can be 
checked directly from ( | C.3| ). 

If we take into account the vanishing of the off diagonal elements of 



ip to be able to replace in (| C.3) the remaining T's by W, the common 



denominator D{z) in ( C.3| ) can be written as: 



D{z) = z'^ - z{Wio,io + Woi,oi)-{uJi-aJof + {uJi-uJo){Woi,oi-Wio.io) 
+ VFio.ioWoi.oi - VFio.oiWoi.io ( C.5) 
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Using ( I 3.1[ - | 3.^ ) the denominator D{z) takes a simpler form: 

D{z) = z'^ - z{Wio.io + Woi.oi) - K - ujof 

+ {uji - ujo)iWoi.oi - VFio.io) ( C.6) 

The analysis of the poles of the diagonal-diagonal elements in Ref. [l^] can 
be repeated here for the elements relative to the diplole moments: D^^{z) 
has the two poles for the same values of z for which the elements Raa.aa{z) ( 



4.1 ), for a = 1 and a = 0, are singular, i.e. for the values z = oji+C,—u)q+5 = 
(5io and z = ojQ + 5 — uJi + Q = 5m. We will note respectively by ^lo and ^oi 
the residue of D~^{z) at those poles, with the obvious property ^lo = Aqi- 
The difference with the unreduced formalism is that those residues are no 
longer given by simple products of ^i, ^o defined from the Green's functions 
associated with the Hamiltonian. The elements of R{z) can now be written 
more explicitly as 



Rn 



Z-LVg + UJa -Waa,aa{z) 

D{z) 



ii-aa.a-a{Z} - ^(^) " ^(^' I ^- 

from which the elements of S(t) can now be computed. We have formally 

Sio.io(i) = e-'^^°Mio(2wi + C-2a;o-^-VFoi.oiK + C-wo + ^)) 
+ e-^^oi* Ai {-Q + 8- t^oi.oi {-^1 - C + ^0 + '^)) , 

Soi.oi(i) = e-^^oiMoi(-2c^i-C + 2u;o + <5-VFio.io(-c^i-C + wo + 5)) 

Soi.io(i) = -e-*^i°*^iot^io.io(^i + C-^o + ^) 

hom{t) = -e-*^«i*AiVFoi.oi(-u;i-C + '^o + 'J) 

- e-'^^''AioWoim{^i + C-^^ + S). ( C.8) 

Those relations are symbolized in the main text as: 

T - -(f'\ — p~*'^10*n' - - -L p~*<^0lt/3 _ _ 

v_ -(f\ — p-iSwt^_ _\p-iSoitf>_ _ f c ql 

where the values of a and /? can be read on the preceding equations ( | C.8| ). 
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A is evaluated from the vacuum- vacuum elements of the operator S(t) 
for the value t = 0. The corresponding elements are: 

-^aa.aa — Olaa.aa ~r Paa.aaj -^aa.aa — Cx.aa.aa ~r Paa.aat \ vviUj 

The elements of the inverse A of the A operator can also be computed. 
This requires the inversion of a two by two matrix. The determinant \Ad\ 
of that matrix is: 

I^dI = (aio.io + /3io.io)(aoi.oi + /3oi.oi) 

- (aoi.io + /?oi.io)(aio.oi +/3io.oi)- ( C.ll) 

We can compute partial contributions from (| C.j )) and obtain 

aio.ioaoi.oi - aio.oiaoi.io = -^^o [2(a;i - a;o)(C - ^) 
- 2(wi - ujo)Wiq,iq{uji + C - wo + 5) + (C - ^)^ 

-(C - i) (Woi.oi(o7i + C - wo + 5) + T^io.io(wi + C - wo + 5))] . 

( C.12) 

We know that 5iq satisfies the equation D{5io) = 0, where D{z) is given by 
( I C.6| ). Therefore, we have 

D{5iq) = SIq- 5io{Wio.io{ho) + WQi,Qi{5iQ))-{uJi-UQf 

+ (a;i-u;o)(VKoi.oi(5io)-l^io.io(5io)) = 0. ( C.13) 

Introducing the explicit value for 5\q {5iq = uji + Q — ujq + 5), this relation ( 
C.13| ) becomes 



(C - ~5f + 2(wi - L^o)(C - ^) - 2{oji - wo)l^io.io(5io) 

+(C - ^)(W^io.io('^io) + Wm.m{5io)) = 0. ( C.14) 



A direct comparison with (| C.12) leads to the relation 



aio.ioooi.oi — oio.oiaoi.io — 0, ( C.15) 

and a similar relation holds for the /3's. 

/?io.io/?oi.oi — /?io.oi/?oi.io = 0. ( C.16) 

The determinant \A£,\ takes the simplest form 

I^dI = aio.io/^oi.oi + aoi.oi/3io.io — aoi.io/?io.oi — oio.oi/^oi.io- ( C.17) 
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A little algebra provides the explicit result 

\Ad\ = AiqAqi [(2lji + Q-2UJQ-5) (-2wi - C + 2c^o + <J) 

+ (C - ^) (-C + ^)- 2(wi - ^o) (VFio.io(-^i - C + ^0 + <5) 
-W^i.oi(wi + C-^o + <^)) 
-(C - 5 - C - (5) (VFio.io(-wi - C + '^o + <5) 
+VFoi.oi(u;i + C-wo + ^))] , 

( C.18) 

In a perturbation expansion, |^d| starts as — 4(a;i — ujqY'. The elements of 

--1 
the inverse A are 

~— 1 1 

\-^ Jaa.aa — ^ \0^aa.aa ~t~ Paa.aa/i 

\M 
D Determination of x 

Denoting by |x| the determinant of the x matrix, we have from the inversion 
of a standard 2x2 matrix 

ix )ii.ii = Xoo.ooi— f) 

\x\ 

(X )00.00 = Xll.lll— T) 

Ixl 

(x" )oo.ii = -Xoo.iij— |, 

(x~"^)ii.oo = -Xii.oorT- ( ^-1) 

Using the explicit value for the elements of 0, we get for the diagonal 
and off diagonal elements of the relation (a = b and a ^ b): 

-u / -l^ _ H^ii.ii(O) 



Xii.ii(x )ii.ii - Xii.oo(x )ii.ii 



VFii.ii(0) + W^oo.oo(0) 



/ -u / -u _ W^oo.oo(O) 

Xoo.iiu Jii.oo — Xoo.oolX Jii.oo 



W^ii.ii(0) + m)o.oo(0) 
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Xii.ii(x )ii.oo — Xii.oo(x )ii.oo 
Xoo.ii(x )ii.ii — Xoo.oo(x )ii.ii 



Woo.oo(O) 



VFii.ii(O) + 1^00.00 (0)^ 
Wii.ii(O) 



(D.2) 



Introducing the form ( | D.l ) for the inverse of x iiito these relations ( | D.2D , 
we get 



Xii.iiXoo.oo - Xii.ooXoo.oo — Ixl 



l^ii.ii(O) 



"Xoo.iiXii.oo + Xoo.ooXii.oo — Ix 
-Xii.iiXii.oo + Xii.ooXii.oo = — Ix 



VFii.ii(0) + 1^00.00(0)' 
W^oo.oo(O) 



Xoo.iiXoo.oo — Xoo.ooXoo.oo 



IXl 



H^ii.ii(O) +1^00.00(0)' 
Vt"oo.oo(0) 
W^ii.ii(0) + M^oo.oo(0)' 
T^ii.ii(O) 
W^ii.ii(O) +1^00.00(0)' 



(D.3) 



Introducing the form (D.l) for the inverse of x in (| 5.8|), we get 



Xoo.oo - Xoo.ii 
Xii.ii + Xoo.ii 



Ixl, 



Xii.ii — Xii.oo — 1x1 
Xoo.oo + Xii.oo = 1- 



(D.4) 



These relations enable to express all elements of x in terms of |x| 



Xii.ii — Xoo.oo 
Xii.oo = Xoo.ii 



^(1 + lxl 
^(1-lxl 



We have indeed: 



Xii.iiXoo.oo - Xii.ooXoo.il = ^(1 + Ixl)^ - ^(1 - Ixl)^ 



(D.5) 



(D.6) 



Introducing those value s (| D-5| ) into ( D-3| ), we get two relations (the third 
and fourth relations in ( | D.3D provide trivially the same relation as the first 
two ones) 



(1 + |X|)'-(1 + |X|)(1-|X|) = 4|x| 
-(1-|X|)' + (1 + |X|)(1-|X|) = 4|x| 



^11.11 (0) 



T^ii.ii(0) + W^oo.oo(0) 

^00.00 (0) 
T^ii.ii(0) + W^oo.oo(0) 



(D.7) 
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Let us add and substract the two relations ( D.7| )- We then get: 

4|x|, 



(i + |x|)^-(i-W)' 
(i + |x|)' + (i-W)'-2(i + |xl)(i-W) 



Mx 



Wn.njO) - Woo.ooiO) 
VFii.ii(0) + VFoo.oo(0)' 

(D.s; 



The first relation is trivially satisfied and we are left with the condition 



4|xP = 4|x 



^^ii.ii(0)-T^oo.oo(0) 



therefore 



IXl 



VFii.ii(0) + W^oo.oo(0)' 
W^ii.ii(0)-VFoo.oo(0) 



(D.9) 



( D.IO) 



W^ii.ii(0) + VFoo.oo(0) 

From that expression, we get the explicit values for the first serie of elements 
of X reported in §5. 

We now turn to the other elements Xaa bb- ^^ ^^^ expression x(exp —i^t)x^^ 
is diagonal and we get 



-iet 



7 y Xaa.bb^ 



KX Jbb.cc- 



( D.ll) 



Yl ^aa.bb ( ^ 
^ \ / bb.cc 



Direct comparison between ( | D.ll ) and ( 4.4 ) provides us the conditions: 

Xaa.loiX )l0.cc = 

Xa-a.OliX-'Wcc = T.^a-a.bb(^"] _ ■ ( D.12) 

(, V / bb.cc 

We introduce the value of the inverse of the operators x ^-rid A. The form 
of the inverse depends on their diagonal or off-diagonal character and the 
two possible cases have to be distinguished 



Xaa.lOXOl.Ol 



Xaa.lOXlO.Ol 



Xaa.OlXOl.lO 



Ixl 



\Ad\ 

\x 



\Ad 

1x1 



l^i^l 



(aaa.io(aoi.oi + /?oi.oi) — aaa.oi(aoi.io + /?oi.io)) ) 

( — aaa.lo(oi0.01 + /^lO.Ol) + "aa. 01 (oiO. 10 + /^lO.lo)) , 
(/?aa.lo(a01.01 + /5oi.Ol) — /3aa.0l(a01.10 + /?01.1o)) > 
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Xaa.OlXlO.lO 



\x\ 



(— /?aa.lo("l0.01 + Ao.Ol) + /3aa.0l(ai0.10 + /^lO.lo)) 



ID I 



( D.13) 
The last two equations provide us the same condition as the first two ones, 



if we take into account the relations ( C.15 - C.16 ). The first two conditions 
( I D-13| ) provide, for a = 1 and a = 0, four relations that are simplified using 
the relations (| ClSfj C.16|) 



Xio.ioXoi.oi 



Xoi. 10X01.01 



Xio.ioXio.oi 



Xoi.ioXio.oi 



IXl 



\Ad\ 

Ixl 
Ix 



\Ad\ 

1x1 

'\Ad\ 



(aio.io/3oi.oi — «io.oi/?oi.io) > 
(ooi.io/^oi.oi — aoi.oi/^oi.io) > 
aio.io/5io.oi + aio.oi/3io.io) > 



(-aoi.io/^io.oi + aoi.oi^io.io) • ( D.14) 



Hermiticity of p^^ (t) implies the obvious property pi^ (t) = /Ogi (t) that we 
impose also on the matrix p^^^ and the following conditions have to hold: 



Xio.io 



Xoi.oi) 



Xoi.io - Xio.oi- 



D.15) 



Therefore, we are looking for a solution to ( D.14) under a form that incor- 
porates these conditions: 

Xio.io = xe"^, Xoi.oi = xe~"^, xio.oi = ye*'^, Xoi.io = ye~''^, ( D.16) 



and we have \x\ = x"^ — y"^. Introducing those values into (| D.14 ), we get: 
2 x2 - 7/2 



1^ 



(aio.io/5oi.oi — aio.oi/?oi.io) 



D\ 



xye 



i{ip-ii) 



x^ - y2 



D\ 



xye 



-i{^-il)) 



1^ 

X^ — 2/2 

' \Ad\ 
x2 — y2 



(«oi.io/5oi.oi — "oi.oi/^oi.io) 



(— aio.ioAo.oi + aio.oi/^io.io) 
(— ooi.ioAo.oi + aoi.oiAo.io) 



D.17) 
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The first and fourth relations ( | D-17| ) are obviously identical when the ex- 
pression of I^dI from ( I C.17| ) is taken into account: their substraction pro- 
vides an identity . The second and the third ones are complex conjugate of 
each other, and provide the value of the phase difference (99 — -0). The first 
relation leads to 

y (aio.io/^oi.oi — aio.oi/3oi.io) = —x (aoi.oi/?io.io — aoi.io/3io.oi) ( D-18) 

and determines y in function of x. A further relation is needed to fix the 
value of X. A remaining indetermination is not new inside the context of 

subdynamics.|14] It does not affect the evolution equations but concerns the 

—P - 

relation between p and p. We can connect, for instance, the value of |x| 

with that of |^d| in a "usual way": |xP = \-^d\ or introduce a criterion that 
ensures the positivity of the density operator. We will not elaborate here 
further on this point. 

From the point of view of a perturbation expansion, the elements x, 
yjOio.io and /?oi.oi take the value 1 for y — > 0, while the other elements 
behave as V'^. 

The relation ( D.18| ) can be made more explicit by replacing the a's and 
the /3's by their formal value. 

X^ [(C -M-S)-iC- 5)t^01.0l(<5oi) - (C - 5)1^10.10 (<5io))] 
= -2/2 [(2a;i + C - 2uJo - 5){-2loi + C + '^^o - S) 

- (2a;i + C - 2wo - ^)l^io.io(5oi) - (-2a;i + C + 2ujo - -5)Woi.oi(5io)] ■ 

( D.19) 



E The presence of an incident field 

E.l A passive incident field 



To evaluate ( | 6.1| ), the matrix elements of the resolvent R{z) are required in 
terms of irreductible operators, starting from the perturbation expansion 



Ra\b.c\d{z) = ^{R{z) 
n=0 ^ 

and a similar expression for Rab\.cd\{z)- 



LvR (z) 



aXb.cXd 



(E.l) 
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The irreductible operators WaXb.cXd and Wabx.cdX for the complete Uou- 
vilUan L are defined through ( 2.19| ): 



WaXb.cXdiz) = Y^ 
n=0 



LvR (z) 



U 



aXb.cXdiirr) 



(E.2) 



WabX.cdx{z) = X! 
n=0 



LvR{z) 



L 



( E.3) 

abX.cdX(irr) 

The irreductibihty condition holds with respect to the vacuum chosen, namely 
the states involving the atom and the incident or emitted photons. 

Since the field line A is not involved in any vertex and plays no role in 
determining the irreductibihty condition, we have the obvious property: 

WaXb.cXd{z) = Wab.cd{z - i^x), WabX.cdx{z) = Wab.cd{z + UJx). ( E.4) 

Those elements enable to write a compact form for the relevant elements of 
the resolvent R: 



RaXb.cXd{z) — Rab.cd{z — OJx), 



R 



abX.cdX 



[z) = Rab.cdiz + UJx) ( E.5) 



The poles that have been considered for the computation of the atomic part 
of S(t) are the poles at z = 0, z = 9, z = 6io, z = Sqi that are present in 
Rab.cd{z). These poles are now shifted in (6.6) by ±0;^ and their residue is 
the same. We use the notations previously introduced: 



-'aXa.bXb 



it) 



e-'^^^'aaaM + e-*(^+"^)*/?. 



(E.6) 



The values of the a's and the /3's can be obtained by identification with 
formulae ( 3.12| ), ( | A3| ) , ( | A8| ), ( | A9| ) or they can be read in the expression 
(["AIsI) for A in §3. 

For the off-diagonal elements, we have 



-'aXa.bXb 



it) 



^ ^aa.bb ^ ^ 



-i{5oi+uJx)tn 



In a similar way, we have 

^aaX.bbxit) 



e^'^^^Oaa.bb + e" 



''aaX.bbX 



it) 



^-i{Sio-^^)t^ 



f-'aa.bb 
-i{Soi-uJx)tf 



_1_ p—'-V"Ul—^\)''Q 
aa.bb ~ ^ r^aa.bb 



(E.7) 

(E.8) 
(E.9) 
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A direct identification is possible: 

-^ab.cd -AabX.cdX = -^ab.cd 

^\Oiaa.bb + {0 + ^\)l3aa.bb 

(5lO + WA)aas.66 + ('^01 + ^x)Paa.bb 
—^XOiaa.bb + {0 — UJx)Paa.bb 



-^aXb.cXd 

(M) 
(el) 
(M) 



( E.IO) 



aXa.bXb 
aXa.bXb 
aaX.bbX 



@A 



--1 



aaX.bbX 



{6iQ - uj\)a^a.bb + ('^01 - '^a)A 



aa.bb 



E.lll 



The inverse A of the operator A is known from the previous section thanks 
to the identification ( | E.IO ). 

We have the obvious property {lab.cd is 1 when a = c and b = d and 
vanishes for the other possibihties) 



(eX) ^^ ^^ = cuxAab.cd+ieA) 

\ J aXb.cXd \ J 

^^~^.bX..dX = -A^--+(^^) 

from which the relations (| 6.2) can be deduced. 



ab.cd 
ab.cd 



( E.12) 



E.2 One absorbed incident field line 

The perturbation expansion of the resolvent R{z) can be written as 



Rab.cXd{z) = ^{R{z) 



n=0 



-0, 



LvR{z) 



( E.13) 



ab.cXd 



A similar expression holds also for Rab.cdx{z) but only one kind of elements 
will be displayed in details. Rab.cXdi^) has to be expressed in terms of the 
irreductible operators Wab.cXd defined in ( | 2.19 ) for the complete liouvillian 

I: 

-0 



Wab.cXd{z) = Y^ 
n=0 



LvR (z) 



U 



( E.14) 



ab.cXd{irr) 

Those elements enable to write a compact form for the relevant elements of 
the resolvent R: 



Rab.cXd{z) — / , / ,Rab.ef{z)Wef.gXh{z)RqXh.cXd{z)- 
e/ gh 



( E.15) 
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Since we have (the hne A is passive through the resolvent and the property 
(I E.5 ) is used to get 



cXd{z) =^Y1 Rab.ef{z)Wef.g\h{z)Rgh.cd{z - i^x) ( E.16) 

e/ gh 

As in §3, in absence of photons, Rab.ef{z) = Rab.ef{z) The elements of the 
resolvent in ( E.16) are therefore the same one's that have been studied 



previously. 

The poles to be considered for the computation of S(t) are the poles at 
z = 0, z = 9, z = 6io, z = 6oi that are present in Rat.efiz) and Rgh.cd{z)- 
These poles are well defined, as we have seen in the preceding sections. 
Moreover, we have seen in Ref. [^] that it is possible to write the elements 
of ^ as a sum of expressions such that each of them contains only one 
relevant pole. The computation of the residue does not therefore present 
conceptual problems. The residues can also be evaluated directly from the 
expression ( |E.16| ). 



^aa.dXdit) = X! Oicia.bbWbb.cXciO)Rcc.ddi-^x) 
b,c 

+ e-''' Yl Paa.bbWbb.cX-c{0)Rcc.ddi^ " '^a) 
b,c 

+ e-^(^io+"^)* Y. ^aa.bbiSio + uJx)Wbb.cXc{Sio + u:x)a,,^ad 

b,c 

+ e-^(^oi+-^)* Y ^aa.bbiSoi + UJx)Wbb.cXciSoi + OJx)Pc-c.dd 
b,c 

( E.17) 

Let us make the following comments for the correct computation of Rcc.ddi~^x), 
Rcc.ddi^ — idx) Raa.bbi^iQ + "^a) and the similar expressions. Displaying the 
two relevant poles and their residue, we can indeed write 

11 

RccAdiz - UJx) = acc.dd- ; T— + Pcc.dd- ; T- + ^ccAdi^ - ^x) 

Z-UJx- Oio Z- UJx- Ooi 

( E.18) 

where the remaining function r^^^^{z — uJx) is regular. When we take the 
residue at the point z = and z = 9, we merely replace z in that expression 
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by its corresponding value. 

KcAdi-^x) = (^ccAd ■ T- + Pcc.dd ■ T— + ^cc.dd(-^A) 

-LOx - oio -^x - ooi 

^cc.ddi^ - ^x) = Oicc.ddE 1 ^ l^cc.dd-S 1 ^ '^cc.ddi^ - ^X) 

( E.19) 

Therefore, we do not have to consider a possible deferred analytically contin- 
uation with respect to the integration variable ti;;^.Q For consistency, when 
taking the residue at z = 5io + wx, we first write Raa.bb{z) as: 

Raa.bh{z) = aaa.bb- + Paa.bb 3 + raa.bb{z) ( E.20) 

z z — a 

and we have: 

Raa.bbiho + l-^x) = Oiaa.bb 



5w + ^X 



+ Paa.bb-, : ^ + raa.bb{5lQ + ^a) ( E.21) 



We have therefore by direct identification from ( E.17| ): 



"^aa.dXd - X! (^<ia.bbWbb.cXc{0)Rcc.ddi-^x) 
b,c 

+ E Paa.bbWbb.cXc{0)R,s,dd{^ " UJx) 
b,c 

+ 51 ^aa.bbiho + ^x)Wbb.cXc{Sio + (^x)acc.dd 
b,c 

+ Y. Raa.bbiSoi + UJx)Wbb.cXc(.Soi + UJx)Pcc.dd ( E.22) 



b,c 



(§1) = ~eY,(3aa.bbWbb.cX-c{0)Rcc.dd{^-^>) 

\ / aa.dXd — 

b,c 

+ ((^lO + ^a) X! ^aa.bbiha + ^x)Wbb.cXc{ho + ^x)oicc.dd 



b,c 

+ (<5oi + Wa) ^ Raa.bbi^Ol + ^x)Wbb.cXc{Soi + UJx)Pcc.dd 
b,c 

( E.23) 
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The inverse A of the operator A can be computed easily. From AA = I, 
we get 

= (aa'') _ = Y.^<^^-bbi^~')bb.dXd + T.^^^-cXci^~')cX-c.dXd ( E.24) 

\ / aa.dXd ^ c 

Since we have obviously (the field line A is purely passive and that case has 
been treated in the previous subsection): 

CA~')cXc.dXd = C^~')cc.dd (E.25) 

and since ^^add' (^ )cc.dd are known from §4 while Aaa.bb, (A )aa.bb are 
known from §3, we have 



(^ )aa.dXd = -Y.C^ )aa.bb\b.cX-cCA )cc.dd ( E.26) 



Those expressions enable the computation of the following elements of 0: 
^aa.bXb^ ^bb.aXa^ ^ aa.bbx^ ^bb.aaX' We have indeed: 

Qaa.bXb = (iQA)A'') 

\ / aa.bXb 

~ /_^{^^)aa.cc{A )^^jjxb + 2_^{QA)aa.cXc{A JcXc.bXb 
c c 

( E.27) 



Using the expressions (| E.24[) and (| E.26D for {A )cXc.bXb and {A )cc.bXbJ 
we get 

^aa.bXb = ^^i&Ajaa.cciA ) cc.ddAdd.eXei^ Jee.bb ~^ Z_^i®^)aa.cXciA )cc.bb 
c,d,e c 

= Y.^-^-ddhd.eXe{A~^)ee.b-b + T.(^^)^--cX-ci^~^)c-c.bb (^-28) 

d,e c 

The last result is reproduced in ( | 6.6| ). 
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